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16090 THIN SECTIONS COMPRESSION 


KEY WORDS: Bulking; Compression; Compression tests; Nonlinear 
programming; Plates (structural members); Struts; Thin plates 


ABSTRACT: The finite strip method developed Cheung extended include the 
nonlinear response imperfect plate strips under axial compression. Local buckling 
and postbuckling phenomena are included the formulation. The convergence and 
accuracy the method tested against accepted solutions for uniformly compressed 
plates with different boundary conditions. Studies are made the nonlinear response 
box and I-section columns. The application the method the interaction local 
and Euler buckling straight columns with imperfect plates presented. 


REFERENCE: Hancock, Gregory J., “Nonlinear Analysis Thin Sections 
Compression,” Journal the Structural Division, ASCE, Vol. 107, No. ST3, Proc. 
Paper 16090, March, 1981, pp. 455-471 


16117 BRIDGE DESIGN AND REGIONAL ESTHETICS 


KEY WORDS: Bridges; Design; Design practices; Esthetic properties; 
Esthetics; History; Regional geography; Structural design 


ABSTRACT: Leading bridge designers have developed individual styles within limited 
locales, and those styles demonstrate similarities personal esthetic ideas. 
examination the works six major bridge designers characterizes these ideas and 
illustrates the striking fact that each case these designers have done nearly all their 
major work in. one well-defined region remarkably small area. The six bridge 
designers are: Thomas Telford (1757-1834), John Roebling (1806-1869), Gustave Eiffel 
(1832-1923), Robert Maillart (1872-1940), Othmar Amman (1879-1965), and Christian 
Menn (born 1927). Their bridges are least good as, and probably better than, 
any their time; there definite connection between the high quality their work 
and its restriction samll well-defined geographic region. Their works illustrate 
series individual styles that spring from local conditions rather than international 
style. 


REFERENCE: Billington, David P., “Bridge Design and Regional Esthetics,” Journal 
the Structural Division, ASCE, Vol. 107, No. ST3, Proc. Paper 16117, March, 
1981, pp. 473-486 


16127 TEST WELDED TUBULAR JOINTS 


KEY WORDS: Bayes theorem; Cycles; Fatigue (materials); Fatigue tests; 
Residual stress; Statistical analysis; Strength (mechanics); Tube joints; 
Welded joints 


ABSTRACT: Low-cycle fatigue tests were performed for various sizes X-type 
welded tubular joint specimens. The relationship between the repeated load and the 
number cycles failure estimated modifying the test data based the 
Bayesean theorem. The size effect and the influence residual stress, which called 
the annealing effect, are also considered. The statistical approach that described can 
applied the fatigue test data which are stongly influenced the scatter 
welding defects. case low-cycle fatigue welded tubular joint specimens, the 
curve can assumed straight line the regular scales. Concerning the size 
effect low-cycle fatigue, the relative strength welded tubular joint specimens 
shows downward trend with increase the size the specimen. 


REFERENCE: Baba, Shunsuke, Arizumi, Yasunori, and Naruoka, Masao, 
Fatigue Test Welded Tubular Joints,” Journal the Structural Division, ASCE, Vol. 
107, No. ST3, Proc. Paper 16127, March, 1981, pp. 487-505 


16128 BUCKLING WEB POSTS PERFORATED BEAMS 


KEY WORDS: Analysis; Beams (supports); Buckling; Elastic properties; 
Girders; Models; Openings; Perforating; Plastic properties; Stability; 
Statistical analysis; Steel beams; Structural analysis; Webs (supports) 


ABSTRACT: The buckling the web post investigated for steel I-beams containing 
adjacent rectangular web openings, and method presented for determining the 
buckling moment. expression derived for the effective buckling length web 
post, and its validity demonstrated means tests scale models. The elastic 
buckling moment the post could then found from elastic buckling theory, and 
used obtain the actual buckling moment from inelastic buckling analysis due 
Nethercot and Trahair. Analysis practial beams shows that web post buckling can 
affect the ultimate flexural strength perforated plate girders, but not significant 
for rolled sections. 


REFERENCE: Dougherty, Brian K., “Buckling Web Posts Perforated Beams,” 
Journal the Structural Division, ASCE, Vol. 107, No. ST3, Proc. Paper 16128, 
March, 1981, pp. 507-519 


16121 BUCKLING COOLING-TOWER SHELLS 


KEY WORDS: Boundary conditions; Buckling; Cooling towers; Shells 
(structural forms); Stability; Structural analysis; Structural engineering; 
Tests; Walls 


ABSTRACT: The equations the linear buckling theory are given tensor notation, 
and the method their solution explained for the case shells revolution under 
axisymmetric loading. Partial differential equations are reduced system ordinary 
differential equations introducing trigonometrical functions for the unknowns the 
problem. For integration, the differential equations are transformed into finite 
difference equations using the Hermitian (multilocal) method. Available test results are 
calculated and good agreement achieved. Parametric studies with different 
boundary conditions demonstrate the sensitivity the buckling loads the statement 
made for the buckling shape regular trigonometric configuration. The so-called 
classical boundary conditions prove very sensitive this respect. 


REFERENCE: Almannai, Ali, Basar, Yavuz, and Mungan, “Aspects 
Buckling Cooling-Tower Shells,” Journal the Structural Division, ASCE, Vol. 
107, No. ST3, Proc. Paper 16121, March, 1981, pp. 521-534 


16142 EDGE STIFFENED CANTILEVER SLABS 


KEY WORDS: Bending; Box girders; Cantilevers; Concentrated loads; 
Design criteria; Elastic analysis; Slabs; Structural analysis 


ABSTRACT: realistic analysis edge stiffened cantilever slabs having linearly 
varying thickness one direction (and subjected concentrated loads) usually 
done rigorous methods that are either computer based are too complex for 
everyday design office use. Semigraphical manual methods are presented which provide 
comparable those given rigorous methods. The methods provide 
longitudinal moments the slab and also the maximum hogging and sagging moments 
the edge beams. The charts are based grillage analogy analyses slabs and 
therefore are subjected the same limitations the grillage analogy method, which 
demonstrated conform closely the classical plate solution for plates varying 
thickness. 


REFERENCE: Bakht, Baidar, “Simplified Analysis Edge Stiffened Cantilever 
Slabs,” Journal the Structural Division, ASCE, Vol. 107, No. ST3, Proc. Paper 
16142, March, 1981, pp. 535-550 
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16140 EXPERIMENTAL EVALUATION COMPOSITE ACTION 


KEY WORDS: Beams; Comparative studies; 
Deflection; Diaphragms (mechanics); Experimental design; Fasteners; Joists; 
Layers; Loads (forces); Lumber; Shear strain; Shear stress; Stress 
(mechanics); Theoretical analysis 


ABSTRACT: Experimental results for layered lumber diaphragm are presented and 
compared with theoretical analysis. The influences joists, boundary members, and 
layer discontinuities deflections and strains are evaluated. The theoretical analysis 
overestimates deflections and strains. Layer discontinuities and boundary members 
have insignificant effects while gluing joists decking results 18% reduction 
deflection. Further tests are recommended reach more conclusive results. 


REFERENCE: Itani, Rafik Y., Morshed, Hossain M., and Hoyle, Robert J., 
“Experimental Evaluation Composite Action,” Journal the Structural Division, 
ASCE, Vol. 107, No. ST3, Proc. Paper 16140, March, 1981, pp. 551-565 


U.S. Conversion Factors 


accordance with the October, 1970 action the ASCE Board Direction, which stated 
that all publications the Society should list all measurements both U.S. Customary and 
(International System) units, the following list contains conversion factors enable readers 
compute the unit values measurements. complete guide the system and its 
use has been published the American Society for Testing and Materials. Copies this 
publication (ASTM E-380) can purchased from ASCE price $3.00 each; orders must 


prepaid. 


All authors Journal papers are being asked prepare their papers this dual-unit format. 
provide preliminary assistance authors, the following list conversion factors and guides 
are recommended the ASCE Committee Metrication. 


convert 


inches (in.) 

feet (ft) 

yards (yd) 

miles (miles) 

square inches (sq in.) 
square feet (sq ft) 

square yards (sq yd) 

square miles (sq miles) 
acres (acre) 

cubic inches (cu in.) 

cubic feet (cu ft) 

cubic yards (cu yd) 

pounds mass 

tons (ton) mass 

pound force 

kilogram force (kgf) 

pounds per square foot (psf) 
pounds per square inch (psi) 
gallons (gal) 

acre-feet (acre-ft) 


Multiply 


millimeters (mm) 
meters (m) 

meters (m) 

kilometers (km) 

square millimeters 
meters 
Square meters 
square kilometers 
hectares (ha) 

cubic millimeters 
cubic meters 

cubic meters 
kilograms (kg) 
kilograms (kg) 

newtons (N) 

newtons (N) 

pascals (Pa) 

kilopascals (kPa) 

liters (L) 

cubic meters 


25.4 
0.305 
0.914 
1.61 
645 
0.093 
0.836 
2.59 
0.405 
16,400 
0.028 
0.765 
0.453 
907 
4.45 
9.81 
47.9 
6.89 
3.79 
1,233 
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NONLINEAR ANALYSIS THIN SECTIONS 
COMPRESSION 


Gregory 


INTRODUCTION 


The interaction between local and Euler buckling thin-walled compression 
member can produce adverse effect its strength. The effect most 
pronounced when the local and Euler buckling loads are approximately equal. 
Its magnitude depends upon the plate imperfections and lack straightness 
the column. Van der Neut (20) and Graves-Smith (5) first demonstrated 
this mode behavior for box columns using mathematical models which 
included the nonlinear membrane stiffness the plates forming the columns. 
The nonlinear stiffness results from the interaction the plate imperfections 
with the local buckling and postbuckling phenomena. 

The determination the interaction between local and Euler buckling 
thin-walled compression member requires knowledge the nonlinear 
membrane stiffness the member. Consequently efficient method required 
for determining the nonlinear stiffness thin-walled section. The finite element 
method has been used successfully Gallagher (7) simulate the nonlinear 
and postbuckling behavior plates. However, the finite strip method developed 
Cheung (3) has the advantage considerably reducing the computation 
involved for uniform prismatic members. The buckling loads and modes 
prismatic members have been investigated the writer (8) using the finite 
strip buckling analysis developed Plank and Wittrick (13). Recently Graves- 
Smith and Sridharan (6) have produced finite strip method for the post-local 
buckling analysis perfect plates and sections. 

this paper the finite strip method extended include the nonlinear 
response imperfect plate strips under longitudinal compression. The displace- 
ment functions used the analysis differ slightly from those Graves-Smith 
and Sridharan. The accuracy the method tested against accepted solutions 
for plates with different boundary conditions. 

Lecturer Civ. Engrg., Univ. Sydney, New South Wales, Australia. 

Note.—Discussion open until August 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication February 26, 
1980. This paper part the Journal the Structural Division, Proceedings the 
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Studies the nonlinear response box and I-sections under axial compression 
are presented. The effective flexural resistances the box and I-sections under 
axial load are evaluated these studies applying small curvature the 
sections and then calculating the resistance from the resultant moment. The 
flexural resistances determined are compared with those Bijlaard and 
Fisher (1), Graves-Smith (5), and the effective section method based the 
Winter effective width formula (4). The flexural resistances can used 
determine the bifurcation load for over-all buckling perfectly straight columns 
with plate imperfections. 


The finite strip buckling analysis folded plate structure has been represented 
(8) 


which [K] linear component the stiffness matrix the folded plate 
assembly; stability matrix which function the longitudinal 
membrane stresses {8} vector nodal line displacements; and 
load factor. Eq. represents eigenvalue problem which can solved 
give the buckling load factor and the corresponding buckling mode {8} once 
the membrane stress distribution specified. 

Graves-Smith and Sridharan extended Eq. the postbuckling range. They 
assumed that the plate assembly was compressed between rigid frictionless platens 
and included nonlinear terms the stiffness matrix. The resulting set nonlinear 
equations represented 


which [X,(8)] nonlinear component the stiffness matrix which 
linear function the nodal line displacements; nonlinear component 
the stiffness matrix which quadratic function the nodal line displacements; 
strain vector caused the longitudinal compression the plate assembly 
between rigid platens; and load vector resulting from the longitudinal 
compressive strains. For the postbuckling analysis, the stability matrix 
formulated function the longitudinal compressive strains rather 
than stresses since the stress distribution changes from initial uniform 
distribution the postbuckling range. 

Eq. has been extended herein include geometric imperfections the 
plate. the initial value {8}, then Eq. becomes 


which represents set nonlinear equations that can solved for {8} 
any value the strain vector The vector must nonzero there 
solution values longitudinal strain less than the compressive 
strain for local buckling. consequence the plate imperfections, the 
equations will produce finite solution equal when the plate unstrained. 
The numerical methods used solve Eq. for displacements are summarized 
later section. 
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Displacement Fields.—In this paper, only one term the Fourier series 
describing the displacement fields along the length strip considered. 
This will produce accurate solution for loads approximately 1.5 times 
the local buckling load when the wavelength this single term the wavelength 
the local buckle the prismatic folded plate assembly. The buckle wavelength 
can determined along with the buckling coefficient solving Eq. set 
out Hancock (8) and Plank and Wittrick (13), charts such those 
produced Bulson (2). 

The rectangular strip shown Fig. assumed compressed and bent 
its plane rigid frictionless platens located its ends. The longitudinal 
compressive strains are assumed nodal lines A,B,, A,B, 


1.—Membrane Displacements Strip: Hookean Displacements; (b) Geomet- 
ric Nonlinear Displacements 


respectively, and vary linearly between these two lines. The resulting 
displacement field the plane the strip (membrane displacements) defined 
These membrane displacements consist two components. They are those 
resulting from Hookean shortening [u,,, Fig. 1(a)] and those geometric 
nonlinear components resulting from plate flexure Fig. 1(b)]. The 


total membrane displacements are simply the sum the two components, that 


The membrane displacements corresponding the Hookean deformations are 
given 


UR y 
L LS 
Parabola 
(a) 
y 
2 
(Sine)? curve 
B> 
(b) 
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the strip; Poisson’s ratio expansion the strip; and the coefficient 
initially. proposed Timoshenko and Gere (17) [sections 9.13, equation 
the term can used allow for the change Poisson’s ratio expansion 
which results from the change membrane stress distribution the nonlinear 
range structural response. 

The geometric nonlinear components membrane displacements which occur 
result the flexural displacements the strip are given 


coefficients solved within the finite strip analysis. The function has 


Sine Curve 


FIG. Displacements Strip 


been chosen sine squared function eliminate shear straining the ends 
the strip. 
The displacement field normal the strip (flexural displacements) consists 


the v-component the y-direction shown Fig. The mathematical 
formulation this displacement given 


which sin This function identical with that used Cheung 
(3). The terms (a,, a,) are the polynomial coefficients solved 
within the finite strip analysis. The component the initial out-of-plane 
imperfection the plate. takes the same form The coefficients 


define the geometrical imperfections the strip and are determined 
prior the analysis. 


y 
8,2 
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The vector nodal line displacements has been chosen exclude the Hookean 
displacements resulting from the longitudinal compression. Its components are 
given 


The displacements are those the central section 
the strip shown Fig. and Fig. However, the displacements 


are those the quarter point section the strip shown, Fig. 


The eight unknown coefficients Eqs. and can represented the 
vector {a} 


The vector {8} can related the vector {a} appropriate substitution 


for (x,z) Eqs. and The resulting set equations can represented 
matrix notation 


{8} 
Solution Eq. for {a} produces 


The matrix [C] given Reference 

Nonlinear Strains.—The nonlinear strain vector, including initial plate imper- 
fections, the same that used Hu, al. (10). includes only those 
nonlinear terms associated with plate flexure. represented 


Vo) 


(=) 
— + 
The linear and nonlinear components the strain vector can derived 
function the nodal line displacements {8} substitution for 
from and into Eq. 12. The resulting linear and nonlinear 
components the strain matrix are given Reference 


Plate Theory.—The plate theory used herein the same that presented 
Cheung (3). This the orthotropic theory derived Timoshenko and 
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Woinowsky-Krieger (18). The stress vector derived using this theory 
M,| 


M,, Gt’ az 12 

E,(€,+v,¢,)/(1 —v,¥,) 


such that The formulation Eq. into property matrix 

Development Stiffness Matrix.—The nonlinear stiffness matrix was developed 
using the theory proposed Powell (14). This theory based upon the virtual 
work principle applied infinitesimal deformations from equilibrium deformed 
position loaded structure. The nonlinear strain terms (Eq. 12) are included 
the calculation the stress and strain when evaluating the internal work. 
The original fixed frame reference has been used for the nodal line displace- 
ments strip when calculating the nonlinear stiffness matrix. The resulting 
set stiffness equations including the geometric imperfections and longitudinai 
compressive strains given Eq. The complete derivation well the 
individual matrices Eq. are set out Reference 

The flexural component the stiffness matrix [K] identical with that 
developed Cheung. The membrane component [K] differs from that 
used Cheung because the use the sine squared function the membrane 
displacement field. The matrix the same the flexural stability 
matrix developed Przemieniecki (15), except that compressive strain rather 
than stress has been used the formulation. 

The [X,(8)] matrix linear function the plate flexural displacements. 
produces the nonlinear coupling between membrane and flexural behavior. 
can broken into two components 


The first component relates the flexural actions the membrane displacements. 
The second component relates the membrane actions the flexural 
displacements. Consequent upon the factor 1/2 applied the second term 
Eq. 14, the resulting stiffness matrix not symmetric. This accordance 
with the theory Powell. The component the stiffness matrix 
quadratic function the plate flexural displacements. results 
modification the flexural stiffness the strip the flexural displacements 
increase. 

Solution Nonlinear Equations.—The simultaneous equations represented 
Eq. are nonlinear and require iterative solution. The Newton-Raphson 
method presented Gallagher (7) for nonlinear finite element problems 


has been used herein. Rearranging Eqs. and imbalance load vector 
produces 


THIN SECTIONS 


nonlinear analysis must proceed estimate and improve upon 
the estimate succeeding iterations. expanding Eq. Taylor series 
and disregarding higher order terms, then after iterations 


6) -[ -{ F,} 


because are seeking {0}. 
The matrix called the tangent stiffness matrix and can 


derived from Eq. 15. this paper, approximate tangent stiffness matrix 
given 


has been used. 

The approximate tangent stiffness matrix symmetric. can derived 
simply from the components [K], which constitute 
the nonlinear stiffness matrix. has been found produce rapid convergence 
the solution the nonlinear Eqs. for loads twice the local buckling 
load. 

Stresses.—The transverse membrane stresses can derived from 
12, and 13. For plate strip which free expand laterally the 
integral over the length the strip zero. The numerical value 
(in can determined equating the integral over the area 
the strip zero. The value determined for each strip after each 
iteration the analysis. 

The longitudinal membrane stress any point strip can also 
derived from Eqs. 12, and 13. The longitudinal stress can integrated 
over the cross-sectional area the strip give the load the strip. The 
values the total load the center and ends folded plate system were 
found agree least significant figures after convergence the solution 


any strain level. This confirms that the solution satisfies longitudinal equili- 
brium. 


The method described the foregoing section has been programmed for 
digital computer. The accuracy and convergence the method has been 
investigated comparing the solutions for square plates with those Yamaki 
(21). His solutions are generally regarded precise values load equal 
least times the local buckling load. 

Two different boundary conditions have been investigated. The first considered 
simply supported square plate whose longitudinal edges are free wave 
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[Yamaki The second condition studied that where the longitudinal edges 
are built-in prevent out-of-plane edge rotations but are free wave in-plane 
Case 

both cases, the difference between the lateral deflections successive 
cycles was less than 0.1% the total deflection after iterations. The same 
degree convergence also applied the longitudinal load carried the plate. 

Simply Supported Plates.—The curve load versus compressive strain calcu- 
lated for the simply supported plate shown Fig. result symmetry 
was necessary analyze only half the width the plate using 
equal width strips. The solution Van der Neut (20), which assumes that 
the postbuckled deflected shape the same the buckling mode, has also 


Van der Neut 
0-408 


@ Van der Neut 
Finite Strip 

© 4 strips over half plate 
6 8 strips over half plate 


Longitudinal edges free wave 


FIG. 3.—Load-Strain Curve Simply Supported Plate 


been plotted Fig. deviates from the Yamaki solution loads above 
1.1 P,. The straight line with slope 0.408 times the initial stiffness (21) represents 
the postbuckling stiffness perfect buckled plate assuming change 
deflected shape from the buckling mode. The Van der Neut solution asymptotes 
this line. The finite strip solution for strips over half the plate agrees 
closely with the Yamaki solution 1.8 P,. The solution for strips lies 
approximately half way between the Van der Neut and Yamaki solutions, and 
sufficiently accurate for practical purposes for loads 1.5 

Plate with Built-In Longitudinal Edges.—The ability the method account 
for boundary conditions other than those simple supports has been investigated 
for built-in longitudinal edges. this case, the plate buckles complete 


ST3 


Yamaki 
Pr 
‘ | 
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wavelength over the full length the plate with buckling coefficient 
7.69 (17) 


The buckling mode sine squared function across the width the plate. 

The geometric imperfection chosen has the same shape the buckling mode. 
The load-strain curve shown Fig. which the finite strip solutions 

are compared with the Yamaki solution. The accuracy the resulting solutions 


Finite Strip Analysis 
© 4strips over half plate 
6 8 strips over half plate 


2m x 2m x 3-2mmplate, v ="3 
Longitudinal edges built-in and 
free to wave.Loaded edges 
simply supported a=0.-1 


3-0 
fle, 


FIG. 4.—Load-Strain Curve Plate With Built-in Longitudinal Edges 


similar those for the simply supported plate, which demonstrates that 
the accuracy the method not significantly affected the boundary 
conditions. 


The extension the analysis from single plate prismatic member 
formed from number plates requires consideration the equilibrium 
nodal line forces and the compatibility nodal line displacements the 
plate junctions. Graves-Smith and Sridharan (6) have devoted considerable 
attention this problem. They found that their choice membrane displacement 


< 9 
+ 4 
05) 
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functions limited the types plate junctions which could analyzed. For 
plate junction applicable, the connected plates must free expand 
result the Poisson’s ratio effect without inducing transverse stresses. 
Applicable and inapplicable junctions are shown Fig. this paper, the 
over-all equilibrium condition (assumed under the heading ‘‘Stresses’’ the 
foregoing)—that the integral the transverse stresses over the area the 
strip zero—produces the same restrictions were discovered Graves-Smith 
and Sridharan. 

comparison the edge flexural displacements one strip and the edge 
membrane displacements orthogonal strip shown Fig. can 
seen that although these displacements match the midpoint and ends 
the strip, there incompatibility over the remaining length the strip. 
This incompatibility likely insignificant for loads since the 
membrane displacements the edge plate are generally order magnitude 
smaller than the maximum flexural displacements adjoining plate. Hence 
the membrane flexural coupling between adjoining plates insignificant. The 
major advantage the assumption Fig. 5(c) limit the total number 


Half wavelength 


of local buckle Flexural 


displacements 
of Strip 1 


Membrane 
displacements 
of Strip 2 


FIG. Junctions: Applicable Junctions; (b) Junctions; (c) 
Relationship Flexural and Membrane Displacements 


degrees freedom line junction with consequent simplification 
programming considerations. The nodal line rotations and longitudinal 
displacements adjoining strips are equated line junction. 

Square Box Column.—An analysis was performed box column which 
had uniform plate thickness and its depth equal its width. quarter 
the column between center lines adjacent plates was studied with each half 
plate divided into four equal width strips. The geometric imperfections each 
plate were taken the same for simply supported plate. The directions 
the imperfections were chosen produce compatibility the plate 
junctions. 

The resulting load-strain curve was identical with that shown Fig. for 
simply supported plate whose longitudinal edges are free wave. This expected 
result confirms that the method can applied structural members. 

The effective axial rigidity the box column was computed from the central 
slope quadratic fitted through three adjacent points the load-strain curve 
Fig. These rigidities are plotted against the load Fig. can seen 
that the geometrical imperfections cause decrease from the initial rigidity 
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the axial load increases. high axial loads, the rigidity drops below 
the value 0.408 calculated Yamaki (21), demonstrating the increasing 
difference between the actual postbuckled shape and the buckling mode. 
Graves-Smith (5) has proposed that the bifurcation load for over-all buckling 
thin-walled column with imperfect plate elements can calculated from 
the Euler formula with the flexural rigidity replaced the effective flexural 
rigidity Bijlaard and Fisher (1) also proposed this approach the post-local 
buckling range thin-walled column composed perfect plate elements. 
Consequently necessary determine the effective flexural rigidity 


Axial rigidity (finite strip) 
Flexural rigidity (finite strip) 


Flexural rigidity (Effective section method) 
Flexural rigidity (Bijlaard and Fisher) 
Flexural rigidity (Graves-Smith) 


FIG. 6.—Effective Rigidities Square Box Column 


column section before the bifurcation load for over-all buckling can calculated. 

The effective flexural rigidity given value the strain vector was 
determined applying first axial compressive strain the box shown 
Fig. and then curvature shown Fig. 7(c). The maximum additional 
strain resulting from this curvature was arbitrarily chosen the 
axial strain. The resulting stress distribution was integrated determine the 
bending moment corresponding the curvature. The effective flexural rigidity 
was determined dividing the moment the curvature. 

The calculations were performed with half the box column subdivided 
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FIG. 7.—Model Box Column Under Compression and Bending: (a) Unstrained; 
(b) Axial Compressive Strain (c) Axial Compressive Strain Bending Curvature (p) 
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FIG. 8.—Effective Rigidities 
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into equal width strips. The resulting flexural rigidity values are plotted 
Fig. and are compared with those produced using the effective section 
method (4), the method Bijlaard and Fisher (1), and the solution Graves-Smith 
(5). For loads greater than 1.1 P,, the finite strip solution and the effective 
width method agree closely with the solution Bijlaard and Fisher. Fow lower 
load values, the flexural rigidity sensitive the assumed values initial 
imperfections. The finite strip solution (with 0.1) close agreement 
with the solution Graves-Smith (with 0.11). recent solution Koiter 
and Van der Neut with 0.05 lies very close the finite strip solution 
with 0.1 and Graves-Smith solution 0.11). Koiter and Van der Neut 
have indicated that their solution conservative. 

analysis was performed I-section for which the flange 


Nodal 
in 


Nodal line 
stress in web 


Buckie stress 


flange 


Mean stress in 


20 flange. 


Stowell 
Analysis Only 


FIG. 9.—Longitudinal Stress Distribution 


width equalled the web depth. The plate thickness was constant throughout 
the section. The plate buckling coefficient Eq. 2.6 for this section 
where the value that for the web [Bulson (2), Fig. 7.12]. The buckle 
half-wavelength 1.5 times the web depth (Bulson, Fig. 1.13). The geometric 
imperfections chosen for the section are shown Fig. The imperfections 
the web the I-section were chosen the same the buckling mode 
the web alone (i.e., sinusoidal). The imperfections the flange were chosen 
the same the buckling mode flange outstand alone (i.e., straight). 
The imperfect web and flange were assumed connected 90°. This results 
imperfection the flange tip 1.56 when the imperfection 
the web center line. The imperfections were assumed have the same 
wavelength the local buckling mode. 


2.0 
Buckle stress 
in web 
Mean stress 
web. 
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analysis was performed for imperfection value 0.1. result 
symmetry, only one flange and half the web were considered the analysis. 
The flange outstands and half the web were each divided into equal width 
strips. The axial rigidity, which was derived using the procedure previously 
described for square box column, plotted Fig. Unlike the axial rigidity 
for the box column, its value drops well below 0.408 times the initial stiffness 
loads above 1.2 P,. 

The longitudinal stress distributions the flange and web values 
equal 1.0, 2.0 are shown Fig. The corresponding load values are 0.912 
and 1.371 P,, respectively. The stress distributions the nodal line and 
buckle center line are shown the upper portion Fig. and the mean 
stress shown the lower portion. The mean longitudinal stress the flange 
and the stress distribution obtained Stowell (16) have also been compared 
Fig. strain value 2.0. These are close agreement despite 
the fact that the analysis Stowell applies perfect isolated flange which 
has web which transfer load. 

The flexural rigidity about the minor axis was also calculated using the procedure 
described for square box column. has been compared Fig. with the 
values produced Bijlaard and Fisher for section without plate imperfections. 
The two curves are parallel for loads greater than 1.1 P,, although the finite 
solution lower than the solution Bijlaard and Fisher. The flexural 
rigidity determined using the effective section method Kalyanaraman, al. 
(11) appears very conservative when compared with the finite strip analysis. 


nonlinear finite strip analysis has been shown produce accurate solutions 
for the deflections and stresses imperfect folded plate assemblies under axial 
compression. Only one series term the longitudinal direction necessary 
achieve accurate solution for loads 1.5 times the local buckling 
load. Four strips over half the width simply supported plate the full 
width flange outstand are sufficient for loads 1.5 times the local 
buckling load when determining the load-strain curve plates compression. 

The effective flexural rigidity imperfect thin-walled box column subjected 
loads greater than the local buckling load was found agree with those 
produced Bijlaard and Fisher (1), and the effective section method 
al. (4). For loads less than the local buckling load, the effective flexural 
rigidity was found agree closely with that Graves-Smith (5). The effective 
flexural rigidity imperfect thin-walled I-column subjected loads greater 
than the local buckling load was found agree more closely with the solution 
produced Bijlaard and Fisher (1) than that produced using the effective 
section method Kalyanaraman, al. (11). The latter solution was very 
conservative when compared with the finite strip analysis. 
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The following symbols are used this paper: 


strip plate width; 
square matrix defined Eq. 10; 
strip; 
axial rigidity undeformed section; 
effective value the nonlinear range; 
flexural rigidity undeformed section; 
effective value the nonlinear range; 
imbalance load vector; 
variable factor which allows for the change Poisson’s 
ratio expansion result change membrane 
stress distribution the nonlinear range structural 
shear modulus; 
stability matrix; 
buckling coefficient Eq. 18; 
[K] linear component stiffness matrix; 
see Eq. 14; 
length strip; 
bending moments about axes, respectively; 
twisting moment; 
longitudinal compressive load plate section; 
value local buckling; 
strip plate thickness; 
values resulting from plate flexure; 
values resulting from Hookean shortening; 
rigid body displacement x-direction; 
load vector resulting from initial strains; 


{a} 
{3} 
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Cartesian coordinate axes; 
x/b; 
sin (kz) and (kz), respectively; 


variable coefficients displacement functions; 
maximum amplitude plate imperfection divided 
plate thickness; 

vector coefficients; 

vector nodal line displacements; 

initial value {8} caused plate imperfections; 
longitudinal compressive strain; 

value local buckling; 

(€, 

longitudinal tensile strains directions, respec- 
tively; 

rotation about axis; 

Poisson’s ratio directions, respectively; 

curvature about axes, respectively; 

twist; 

longitudinal compressive stress; 

average value 

value local buckling; 

longitudinal tensile stresses directions, respec- 
tively; and 

shear stress plane strip. 
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BRIDGE DESIGN AND REGIONAL 


INTRODUCTION 


The purpose this paper show how leading bridge designers have 
developed individual styles within limited locales and how those styles demon- 
strate similarities personal esthetic ideas. examination the works 
six major bridge designers will characterize these ideas and illustrate the striking 
fact that each case, these designers have done nearly all their major work 
one well-defined region remarkably small area. 

the thesis this study that the bridges these six men: Thomas Telford 
(1757-1834), John Roebling (1806-1869), Gustave (1832-1923), Robert 
Maillart (1872-1940), Othmar Ammann (1879-1965), and Christian Menn (born 
1927) are least good as, and probably better than, any their time, 
and that there definite connection between the high quality their work 
and its restriction small well-defined geographic region. short, their 
results illustrate series individual styles, rather than international styles, 
that spring from local conditions. Each these styles shows the vision 
single designer strove meet local conditions, and express the 
esthetic potential bridge forms saw it. 

This study summarizes briefly the career each designer, the local conditions 
within which worked, and the individual style each. What sets these 
six designers apart that they all strove consciously design bridges that 
met local physical conditions with little materials possible, that satisfied 
the political utility with little public money possible, and 
that improved the environment visually expressing clearly possible 
their idea appropriate structural form. Thus, this study, the bridges defined 
the best are those which combine, the highest degree, efficiency use 
materials, economy money, and handsome appearance. 

The basis for choosing only six out perhaps the best-known modern 

the April 14-18, 1980 ASCE Convention and Exposition, held Portland, 
Oreg. 


Civ. Engrg., Princeton Univ., School Sci., Princeton, N.J. 
08540. 

Note.—Discussion open until August 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication July 1980. 
This paper part the Journal the Structural Division, Proceedings the American 
Society Civil Engineers, Vol. 107, No. ST3, March, 1981. ISSN 0044- 
8001 /81 /0003-0473 $01.00. 


473 


474 MARCH 1981 ST3 


bridge designers that these six show more clearly than their contemporaries, 
distinct and consistent individual style visual form. short, the choice 
basically esthetic. further true that each the six has probably been 
studied more thoroughly than any his contemporaries, expecially from 
combined esthetic-technical point view. is, however, entirely possible that 
with more detailed study other designers will seen possessing style 
quality equal that any these six studied here. fuller list designers 
the first rank would surely include Brunel (1806-1859), Eugene Freyssinet 
David Steinman (1885-1960), and Ulrich Finsterwalder (born 
1897). There are number other designers still active 1980 who might 
well merit detailed study, well.] 

Although these six designers, except for Menn, are well-known, brief review 
each useful order draw attention common esthetic ideas. Therefore, 
after briefly reviewing the career, locale, and style each designer, shall 
explore their similarities and conclude with few tentative generalizations. 


Telford was born Glendinning, Scotland August 1757, and helped 
build bridge Langhold, three-span masonry arch bridge completed 
about 1778. The first bridge built his own design, also multispan stone 
arch bridge, was completed 1792. His first iron bridge called Buildwas, built 
1795, mile upstream from the 1779 Iron Bridge, had span 130 ft, 
which, though 30% greater than the Iron Bridge span, required only half the 


weight the earlier work. 


FIG. 1.—Craigellachie Bridge 
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His work with metal continued with the Aqueduct built between 
1795 and 1805. This was the most impressive canal structure its day and 
remains today remarkable work. Metal arches, resting masonry 
piers over 100 above the valley, carry the canal over the River Dee. The 
original iron work good serviceable condition today after about 1-3/4 
centuries use. 

Telford’s most modern-looking metal arch bridges are those Craigellachie 
Scotland (150 span, completed 1815); Bettwys-y-coed, called the 
Waterloo Bridge (finished 1815); near Tewkesbury, called the Mythe Bridge 
(170 span, finished 1826); Holt Fleet over the Severn (completed 
1827); and Galton (150 span, built 1829). These are all characterized 
flat lattice cast-iron arches connected the nearly flat roadways above 
thin diagonals. Over about the central third the spans the arches carry 
the roadway directly (1). 

Telford’s two greatest bridges are the 327-ft span Conway Bridge and the 
579-ft span Menai Straits Bridge, both Wales, and both completed 1826. 
They are eyebar chain suspension bridges with light decks and masonry towers. 
Menai was the longest spanning bridge the world until the 1834 completion 
the Sarine Valley Bridge Fribourg, Switzerland. The Conway Bridge 
still service, and has stood well for over 150 yr. The Menai, the other 
hand, suffered much damage from wind and has been rebuilt several times, 
but also still service today. 

For our study, here the flat-iron lattice arches (see Fig. that constitute 
Telford’s mature style medium-span bridges. With the exception Craigella- 
chie, they were built the west England and Wales hilly country cut 
narrow streams. was the region where the Industrial Revolution 
said have begun just the time Telford began practice (2). Telford’s 
arches characterize that revolution their material, i.e., cast iron, and their 
form, i.e., the lattice arch, which would the same basic form for the longest 
spanning steel arches the twentieth century. 


Joun 


John Roebling was born Muhlhausen, Saxony, Germany June 12, 
1806. studied civil engineering the Polytechnic Institute Berlin, finishing 
1828, after which worked the next building roads for the Prussian 
government. May 1831, left Germany with group others for America, 
where they bought land western Pennsylvania and founded Saxonburg 
farming community German immigrants. 

1837, Roebling became American citizen. Bored with farming, took 
job engineer for the state Pennsylvania, building dams and locks, 
and surveying line for prospective railroad route. soon became principal 
assistant the chief engineer the state, and Johnstown, Pa., became 
familiar with the newly-constructed Portage Railroad where long canal boats 
were hauled mountains hemp ropes. Roebling successfully substituted 
iron wire rope for the hemp and established factory for wire rope Saxonburg 
the summer 1841. 

1844, won the competition build his first suspension bridge, carrying 
canal over the Alleghany River; and 1849, when moved his factory 
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Trenton, N.J., was already success, both factory production and 
bridge building. Roebling’s first suspension bridge for roadway was built 
Pittsburgh, over the Monogahela River 1845, and his next major works were 
the 815-ft span Niagara Falls Rail and Road Bridge, completed 1855, and 
the Cincinnati Bridge begun 1856 which was interrupted the Civil War, 
but afterward completed 1866. 

March 1857, Roebling wrote letter Horace Greeley, published the 
Tribune, announcing his intention build East River Bridge. Greeley had 
proposed bridge the Tribune early 1849. Then April 1967, charter 
was granted the New York Legislature Albany, and that September, 
Roebling first presented his plan Brooklyn. His written proposal had begun 
with perhaps his most noted statement which included 


great work art, and successful specimen advanced bridge 
engineering, this structure will forever testify the energy, enterprise, 
and wealth that community which shall secure its erection. 


February 1869, Roebling presented his plans consulting board including 
William Jarvis McAlpine, president the newly-reformed ASCE, and also 
presented his plans Henry Latrobe, son and namesake the architect chosen 
Thomas Jefferson rebuild Washington, D.C. after its burning the British 
1812. June 28, 1869, Roebling’s foot was crushed ferryboat while 
surveying for the bridge and died lockjaw July 22. His eldest son, 


Washington Roebling, became chief engineer for the bridge. Flags Brooklyn 
were flown half mast for Roebling. 


FIG. Bridge 


| 
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The Brooklyn Bridge, under Washington Roebling’s direction, began with 
ground breaking for the Brooklyn Tower foundation January 1870, and 
was Officially opened May 24, 1883. Built essentially designed John 
Roebling, quickly became the landmark predicted, and was officially made 
National Historic Landmark 1964 (3). 

Roebling’s mature style consisted masonry towers, thin deck supported 
vertical suspenders, and cable stays radiating out from the towers (see Fig. 
2). Nearly all his bridges were built the states Pennsylvania, New 
York, and Ohio, and his three major works are all over large open waterways. 


Gustave 


Eiffel was born Dijon, France, December 15, 1832, and after study 
the local schools tried for entrance into the Ecole Polytechnique, but 
failed the examination. enrolled instead the Paris Ecole Centrale des 
Arts Manufactures from which graduated 1855, with degree chemical 
engineering. had intended enter his uncle’s vinegar business, but family 
disputes prevented that, joined construction firm, quickly rising 
chief engineer bridges (4). the 1860s decided establish his own 
firm and 1869 had designed and built series viaducts France’s 
Massif Central that established his mastery metal bridge form. 

1878, Eiffel was regarded the leading builder metal structures 
France, and many, the leader the world. His works included major 
bridges many parts the world, the most important which are: the 1864 
viaduct over the Creuse River Busseau, the four viaducts between Gannat 


FIG. 3.—Viaduct Rouzat 
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THE GARABIT VIADUCT OVER THE TRUYERE.CENTRAL FRANCE. 


FIG. 4.—Garabit Viaduct 


and Commentry 1867-1869 (at Bouble, Bellon, Rouzat, and Neuvial), 
the 1877 Maria Pie viaduct over the Douro Portugal, and the Garabit Viaduct 
the Massif Central, completed 1884 the longest spanning arch the 
world, with 165 (5). 

Eiffel’s designs the Massif Central reflect the problems heavy wind 
loads very high viaducts and his style was widen his towers and arches 
the lateral direction, giving the strong visual impression stability (see Figs. 
and 4). 


Rosert 


Robert Maillart was born February 1872 Berne, Switzerland, where 
his father, Belgian citizen, was banker. studied civil engineering 
the Federal Institute Technology Zurich, and graduated 1894. Ironically, 
one his lowest grades was bridge design, even though regarded 
today one the half-dozen greatest bridge designers since the industrial 
revolution. 

For following his graduation, worked with different civil engineering 
organizations. 1902, founded his own firm for design and construction; 
thereafter his business grew rapidly and expanded far Russia and Spain. 
the summer 1914, took his wife and three children Russia. The 
World War made difficult return Switzerland, Maillart stayed and 
worked Russia until 1918, when his business was liquidated the Revolution. 
Forced flee, returned Switzerland penniless and lonely, his wife having 
died Russia. 

Because these misfortunes, Maillart felt unable take the construction 
business again, and from then concentrated design alone. opened 
office Geneva 1919, adding branches Berne and Zurich 1924. 
During the 1920s, began develop and modify his ideas bridge design; 
until his death 1940, produced over bridge designs extraordinary 
originality. Unfortunately, Swiss municipality would accept his designs for 
prominent urban locations, and most these designs were built rural areas, 
many the mountainous countryside. 

1936, was elected honorary member the Royal Institute British 
Architects although had never officially acted architect any project. 
The 1941 first edition Space, Time and Architecture art historian Sigfried 
Giedion, introduced Maillart wide public the United States. Finally, 
Max Bill’s 1949 book, Robert Maillart, with its photographs and commentary 
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nearly all Maillart’s major bridges, powerfully presented him artist 
the first rank. 

Maillart developed two major new bridge forms reinforced concrete: the 
hollow-box and the deck-stiffened polygonal arch. His early hollow boxes were 
all three-hinged arches, but later used the form straight spans well. 
The 1930 Salginatobel Bridge (see Fig. shows these three-hinged forms, while 
the 1933 Schwandbach Bridge (see Fig. shows the deck-stiffened types. These 
forms were built often narrow valleys the hills and mountains 
Switzerland, nearly all only two the Swiss cantons: the Graubiinden and 
Bern (6). 

Maillart’s style shows extraordinarily-thin polygonal arches vanishingly-thin 
hinged points the hollow-box forms. The concrete fully exposed, and 
the forms provide strong contrast the mountain environment. Most his 
major works got built because they won cost competitions, which nearly 
always proved his designs less costly than alternate designs. The long 
tradition Swiss design/construction competitions allowed Maillart win 
designs that probably would not have gotten had design commissions been 
awarded directly public officials consulting engineers. 


Ammann was born March 26, 1879, Feuerthalen, Switzerland, attended 
schools near and Zurich, and entered the Federal Polytechnical Institute 
1897, graduating 1902. had the same professor for bridge engineering, 
Wilhelm Ritter (1847-1906), did Robert Maillart (7). 1904, Ammann came 
the United States where worked for several different firms ending 
with Gustav Lindenthal from 1912-1923. There worked the 1915 Hellgate 
Bridge, and later Lindenthal’s proposals for Hudson River suspension 
bridge. 

Ammann left Lindenthal 1923, submitted his own design for the Hudson 
River Bridge (now called the George Washington Bridge) late that same year, 
and 1924 was named chief bridge engineer for the Port Authority New 
York and New Jersey, where remained until World War II. After the war 
formed partnership with Charles Whitney. 

Ammann’s major designs were: the 1931 George Washington Bridge with 
span 3,500 ft, the 1931 Bayonne Bridge with arch span 1,652 ft, 
the 1936 Triborough (suspension) Bridge with span 1,380 ft, the 1939 
Bronx-Whitestone (suspension) Bridge 2,300-ft span, the 1,850-ft span 
Throgs-Neck (suspension) Bridge 1961, and the 4,260-ft span Verrazano 
(suspension) Bridge, opened November 1964, just months before Ammann’s 
death September 22, 1965. The remarkable fact that with few exceptions 
(8), Ammann’s long span designs are all New York City and possible 
clear day see them all from the Port Authority’s World Trade Center! 
(He also designed the 1961 Walt Whitman Bridge Philadelphia, well 
many smaller works.) Not any the previous four designers had such narrow 
region which was concentrated his major works. 

Ammann’s style suspension bridges consisted strong vertical towers 
with single cross member the top (see Fig. and light horizontal decks. 


BRIDGE DESIGN 


FIG. 7.—Bronx-Whitestone Bridge 


The only exception this visual pattern the towers the Triborough 
Bridge, but there Ammann did not make the original design, but only modified, 
for the purpose economy, earlier one made others (9). 


Menn 


Christian Menn was born 1927 Meiringen, the Berneroberland, the 
third five children. Because his father’s work, lived various parts 
Switzerland, but his family was from the Graubiinden, the largest and most 
sparsely settled Swiss canton. When Menn was old, went live 
with his uncle farm near the small Graubiinden town Sufers the 
Hinterrhein, just below the San Bernardino pass. began his studies 1939 
the canton school Chur, capital the Graubiinden. 1946, entered 
the Federal Technical Institute Zurich, graduating with degree civil 
engineering 1950, after which entered the military service. 1951, 
contracted tuberculosis and spent months hospital Arosa, where 
began theoretical studies and his doctoral work. After his recovery, went 
back his university and became assistant Professor Pierre Lardy, from 
whom had first learned detail about the work Robert Maillart. 

After receiving his doctorate structural engineering 1956, Menn worked 
Paris the United Nations Educational, Scientific, and Cultural Organization 
buildings designed Pier Luigi Nervi. returned Switzerland 1957, 
and June that year opened his own office Chur. the years that 
followed, designed numerous bridges and buildings, mostly the Graubiinden, 
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and won first prize for several bridge competitions other cantons 
Switzerland. 

1971, Menn was made professor structural engineering his old institute 
Zurich and gave the office Chur. Since then, has been consultant 
most the large bridges Switzerland, and since 1977 has been president 
the commission for the revision the Swiss building code for reinforced 
and prestressed concrete. addition, working with the design offices others, 
had made series major bridge designs. Menn’s Felsenau bridge, built 
1974, now the longest span bridge Switzerland, but will surpassed 
his bridge for the Ganter Valley, completed 1980. 

Menn was strongly influenced the works Robert Maillart and his earliest 
bridges are visually almost copies Maillart’s two major forms (10). However, 
Menn’s style evolved slowly increasing labor costs made arch scaffolding 
and formwork more and more costly. His deck-stiffened arch bridges the 


FIG. and Cascella Bridges 


1960s show extreme lightness and increase spacing between cross 
walls (see Fig. 8), both reduce formwork and take advantage prestressing 
the deck. 

1970, competition was held for very large viaduct built just 
north Bern over the Aare River Felsenau. Seven groups were invited 
submit designs. Menn, head one group, was expected produce 
fine arch design. Instead, because scaffolding for very high works was becoming 
increasingly costly, submitted the plan for fully prestressed cantilever bridge. 
Ulrich Finsterwalder Germany had pioneered this type bridge, which 
built without scaffolding. The forms for small sections are moved out horizontally 
from already built columns, and the sections are cast and allowed harden; 
they are then prestressed tendons, which tie the newly-formed sections back 


BRIDGE DESIGN 483 


into those previously cast. The bridge thus built out itself, using only 
small amount forming and scaffolding from the ground. 

Menn took Finsterwalder’s method and modified the design two significant 
ways. First, used two slender columns each major support, reducing material, 
and giving the visual impression lightness; second, used, for relatively 
wide bridge, only single box with very long slabs cantilevered laterally 
complete the roadway. Unlike most previous long-span cantilevered bridges, 
the Felsenau curved plan, making more complicated (see Fig. 9). 


FIG. 9.—Felsenau Bridge 


Because the economy and elegance Menn’s solution, won first prize, 
and the bridge has since been completed according his plans. The main 
156-m spans, now the longest Switzerland, together with the 60-m columns 
make this bridge impressive sight curves out the woods and across 
the river, high above the ordered suburban communities below. Thus, the 
1970s, the potentials prestressing were being exploited more fully, and arches 
were longer competitive for shorter spans. 


> 
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Even more impressive than Felsenau the design proposed Menn 
1975 for the Ganter Bridge the Simplon Road. Here the canton Valais 
had intended build tunnel. Menn, however, estimated that his bridge would 
cost about half the amount needed for the tunnel, and bids the Ganter 
Bridge proved him right. Completed 1980, this immense work crosses wide, 
deep valley with central span 174 long; the highest the supporting columns 
130 More remarkable than the sheer size the form, prestressed cantilever 
girder which the main prestressing tendons rise well above the girder 


FIG. 10.—Ganter Bridge 


the columns. These tendons are concreted into walls over the central part 
the span give unique profile that fully exploits prestressing (see Fig. 10). 


Each these six designers did evolve personal style which sets their work 
apart from others. addition, that style developed within one relatively small 
region which each designer knew intimately. every case these designers have 
developed new forms and some vastly increased scales. They went beyond 
the state-of-the-art, leading the way new state and new set art forms. 
nearly all cases their designs were publicly debated competition with designs 
others. Each gained prominence making designs which were less costly 
than others. 


Except for Telford, whose early career predates formal technical education, 
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each was well trained the best engineering schools. Each had either direct 
experience builder (Roebling, Eiffel, and Maillart), studied building practice 
firsthand integrate construction ideas directly into design. 

Finally, each designer wrote about esthetics, and each considered central 
design concern. Except for Ammann, none these designers worked together 
with architects their mature works; and Ammann clearly determined the 
form his bridges and did not let anyone else decide major issues appearance 
(11). Thus, each designer believed that the appearance was the full responsibility 
the engineer and that major goal design was handsome form. 

Embury, the author Ref. 11, clearly describes Ammann’s control the 
design, and several other articles details his own ideas which are characterized 
the following: 


belief that more good engineering projects are spoiled false 
deference assumed esthetic considerations than any economic 
requirements the structure whose function partly concealed and 


partly distorted for supposed esthetic reasons will inevitably unsuccess- 
ful. 


The common ideas these six designers can tentatively identified generally 
lightness thinness, exposed and undistorted structure, and form that 
fits the local environment. Each these general ideas needs careful visual 
analysis that considers the works each designer separately. Such analyses 
are beyond this present scope, can only conclude brief clarification 
each idea reference some the figures. 

Lightness thinness means that visually, the new materials industrialized 
metal, reinforced and prestressed concrete, allow the use very thin and open 
sections, and that these designers wanted express that visually (see Figs. 
and 9). This represents denial the centuries-long tradition heavy 
solid closed forms typical most masonry bridges. 

Exposed and undistorted structure means that visually, the new forms were 
the load-carrying parts without any cladding whatsoever. Further, those exposed 
forms were shaped give both greater efficiency and better looking result 
the same time (see Figs. and 10). 

Finally, form that fits the environment means one that fits the foundation 
conditions, the wind forces, the approaches, and the weathering. Fitting the 
environment does not mean making stone facades because rocky terrain, 
any other such superficial visual expression unrelated the structural ideas. 
Thus, the Rouzat towers (Fig. spread reflect visually the overturning 
the wind, the Brooklyn Bridge diagonal stays show Roebling’s concern with 
the vertical effect wind (Fig. 2), and the curved sweep the Felsenau 
(Fig. and the Schwandbach (Fig. 6), fit with the turn the highway. The 
Salginatobel’s flat three-hinged form the steep valley walls and the 
dangerous vertical movements one them (Fig. 5). 

seems that for each these six designers the expression bridge form 
came from the designers: (1) Being restricted well-defined region small 
area and relatively uniform topology; (2) working political context where 
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alternative designs were publicly debated, and often design competitions 
design-construction competitions held for bridges; (3) slowly gaining experience 
through completed works that their best works invariably came later life; 
and (4) gradually evolving personal style that visually clear retrospect. 

All these factors help support the conclusion that these six, along with 
the best their contemporary designers, were and are artists new art 
form, Structural Art, made possible the new materials the industrial 
revolution. Thus, there are three sorts artists the modern world working 
three-dimensional forms: engineer, architect, and sculptor. Each indepen- 
dent artist type, although each can learn from the others and although sometimes 
they collaborate works. review bridge design reveals that structural 
art has developed independently from the arts architecture and sculpture, 
and such, this engineering art holds great promise for the future. 
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FATIGUE TEST WELDED 
TUBULAR JOINTS 


Shunsuke Yasunori and Masao 


Welded tubular truss structures are currently widely employed the field 
marine structures. important problem such structures evaluation 
the strength welded joints under strong repeated loadings the plastic 
range. 

Static and high-cycle fatigue tests welded tubular joint specimens have 
often been reported: for example, Akiyama (1) for T-type, X-type, and K-type 
joint specimens; Kanaya (4) and Toprac (9) for T-type joint specimens; 
and Maeda (6) and Washio (10) for X-type and K-type joint specimens. 
However, far the writers know, report has been published concerning 
low-cycle fatigue tests welded tubular joint specimens. 

The structure will vibrate cyclically the plastic range under the action 
strong earthquake motion strong wave motion, which can regarded 
kind high-stress, low-cycle fatigue process. Experimental and theoretical 
studies the low-cycle fatigue phenomenon have been reported many 
researchers, which the low-cycle fatigue defined fatigue phenomenon 
whose number cycles failure less than 10° 10*. However, fatigue 
tests which less than are extremely scarce even for the simplest 
specimen, such smooth rod. 

this paper, static tensile tests and low-cycle fatigue pull-pull tests 
600) using X-type welded tubular joint specimens are described. Fifty-six 
specimens the same size and the same material, whose height was 1.2 
in.), were prepared both for the static and low-cycle fatigue tests. 

Since welded tubular joint specimens are complicated structures compared 
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with smooth rod, the numbers cycles failure are likely show 
strong scatter. The scatter will probably result from welding defects such 
lack penetration. But the case X-type welded tubular joint specimens, 
the weld lines are long and complicated. Welding defects are therefore hard 
measure along all the weld seams; and even they could measured, 
would impossible evaluate the effects welding defects the fatigue 
strength. Thus, the scatter welding defects not practical factor explaining 
the scatter The writers try substitute the scatter static collapse 
load for the scatter welding defects. 

The writers believe that the large scatter due the scatter welding 
defects (static collapse load) should eliminated considering the fatigue 
property the specimens, and that the relation between the repeated load 
and the number cycles failure which called the curve, should 
recalculated for the particular state which the static collapse load 
all specimens made equivalent applying the Bayesean theorem and the 
concept the weighted mean. 

The effect size the welded tubular joint and the influence residual 
stress, which called the annealing effect, are examined preparing four 
kinds specimens different sizes, which half are annealed order 
eliminate the residual stress. These two effects and the relationship between 
them are considered herein. 


The static tensile test and the pull-pull fatigue test were conducted using 
the X-type welded tubular joint specimens shown Fig. 
Four kinds specimens different sizes—they were similar shape—were 


[| tension 
<> 


FIG. 
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prepared, and are described through case-4 specimens. Let Case-3 
specimen standard-size specimen, which chord diameter 139.8 mm, 
chord thickness 4.5mm, chord length 700mm, brace diameter 101.6 
mm, brace thickness 3.2 mm, and brace length 1200 mm. The dimensions 
the specimens are listed Table which shows that the Case-1 specimen 
had quarter much cross-sectional chord area the Case-3 specimen, the 


TABLE 1.—Sizes Specimens 


TABLE 2.—Numbers Specimens Used Test 


Fatigue 
(3) 


(5) 


Note: Numbers parentheses refer total number specimens, including misses. 


Case-2 specimen had half much chord area the Case-3 specimen, and 
the Case-4 specimen had one and half times much chord area the Case-3 
specimen. The ratios chord areas Case-1, Case-2, and Case-4 specimens 
the Case-3 specimen were not exactly 1/4, 1/2, and 1.5, respectively, and 
also the ratio chord diameter chord thickness (d,/t,) was not strictly 


age sec- 
tional Brace 
Diam- thick- area Diam- Thick- 
Length ter ness A, Length | eter n Among 
9g e g ess 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
300 60.5 421 0.215 500 42.7 0.706 24.2 
114.3 1,218 0.591 900 76.3 0.668 31.8 
700 139.8 1.000 1,200 101.6 0.727 28.0 
950 190.7 3,087 1.431 1,650 139.8 0.733 36.7 
q 
Case 
(1) (6) 
(5) (4) (11) 
(10) (41) (4) (56) 
(1) (5) (4) (11) 
(69) (20) (89) 
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constant for each case specimen, because commercially available standard-size 
tubes were employed. Case-1, Case-2, and Case-4 specimens were used consider 
the size effect comparing them with the Case-3 specimen. 

Four kinds annealed specimens, which were the same sizes the 
through Case-4 specimens, were prepared order consider the influence 
residual stress so-called annealing effect. The numbers specimens 
various sizes are shown Table which the annealed and the nonannealed 
specimens are indicated separately. 

The specimens were made steel tube, whose tensile strength was 
4.02 MPa (41 kg/mm’), that is, whose tensile collapse load was almost 402 
kN, because the cross-sectional brace area was almost 1000 mm’. The welding 
operation was done skilled workman, and the surface the welds was 
not smoothed grinder can seen the practical structures. 

The repeated load was applied constant-load method. was pull-pull 
test, and the minimum pull load was one-tenth the maximum pull load 


TABLE 3.—Tensile Strengths Case-3 Specimens 


kilonewtons 


(1) (2) (3) 


“Mean 5.00. 
290; standard deviation 27.2. 


minimum 9.8 (1000 kgf). The frequency the repeated load was cycles 
per min, which the shortest period for the machine generate exact 
sine-wave loading. 

Details the tests are follows: 


Static tensile tests for Case-3 specimens: Ten specimens were prepared 
for the static tensile test order obtain good estimate the mean and 
the variance static tensile collapse load S,. The results are shown Table 
which the other data coming from the fatigue test are added the specimen 
linear relation between the static collapse load and the cyclic collapse 
load has been reported (3,8,11). The collapse loads are modified related 
chord thickness (see next section), and the modified collapse loads are 
expressed 


4.99 249 250 
4.99 255 256 
5.03 267 265 
5.01 275 
4.97 275 278 
4.97 292 295 
4.99 301 302 
5.03 314 311 
4.97 314 317 
4.97 320 323 
4.97 324 327 
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Pull-pull fatigue tests for Case-3 specimens: Forty-one specimens were 
prepared and divided into five categories according the magnitude repeated 
loads, and successful data were obtained for specimens. Five categories 
repeated loads were follows (numbers specimens are shown 
parentheses): 255 (26,000 kgf) specimens); 245 (25,000 kgf) 
(20,000 kgf) (4). The numbers cycles failure are shown Table 
which the repeated loads are also modified related chord thickness 
well the static collapse load, and are expressed The relationships 


TABLE 4.—Applied Repeated Loads and Numbers Cycles Failure Case-3 
Specimens 


255 (26,000 kgf) 


256 255 253 253 

0.882 0.878 0.872 0.882 0.878 0.872 0.892 
172 
162 184 


244 

0.841 0.828 0.851 
132 227 
199 241 


185 
247 


215 

0.736 0.740 0.747 
410 448 
412 430 


between the nondimensional load and the numbers cycles failure 
are shown Fig. which circular symbols are employed. Fig. 
nondimensional repeated load used ordinate instead 
the stress. Stress may reasonable parameter ordinate the specimen 
has the same cross section along its length; but the case X-type welded 
tubular joint specimen, the stress not constant except part the brace 
distant from the joint. The axial stress, which obtained dividing the collapse 
load whole structure nominal sectional area brace, not easily 


0.872 0.892 
235 (24,000 kgf) 
215 240 233 271 280 332 350 360 
216 (22,000 kgf) 
287 289 
196 (20,000 kgf) 
196 196 196 196 
0.676 0.676 0.676 0.676 
403 473 523 544 
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500 600 


500 600 


FIG. 2.—S-N Relationship all Specimens: (a) Before Modification; (b) After Modi- 
fication 
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TABLE 5.—Fatigue Test Case-1 Specimens 


Note: indicates static collapse loads. 


TABLE 6.—Fatigue Test Case-2 Specimens 


Note: indicates static collapse loads. 


TABLE 7.—Fatigue Test Case-3 Specimens 


Variable 
(1) 
kilonewtons 
N* 


Note: indicates static collapse loads. 


493 
Variable Data 
1.0 1.0 0.976 0.929 0.895 
280 383 153 
1.023 0.919 0.723 0.686 
136 508 441 
Variable Data 
Normal: kilonewtons 167 157 147 128 
kilonewtons 189* 162 156 156 129 
1.017 0.899 0.861 0.834 0.699 
601 160 735 
154 284 217 605 
kilonewtons |211* 169 150 
1.081 0.855 0.760 
133 
215 234 
Data 
332* 292 235 
294 237 
0.990 0.865 0.696 
661 
211 576 
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regarded representative the total stress state structure. Fig. 
both abscissa and ordinate are plotted regular scales instead 
log scales. There fixed rule using log scales for low-cycle fatigue 
problem, and regular scales are the best fitting scales for the test data. 

Static tensile tests and pull-pull fatigue tests for Case-1, Case-2, and Case-4 
specimens: One specimen was prepared for static test, and all five specimens 
were prepared for every fatigue test these specimens. The results are shown 

Static tensile tests and pull-pull fatigue tests for annealed specimens: One 
specimen, again, for static test and all five specimens for every fatigue test 


TABLE 8.—Fatigue Test Case-4 Specimens 


Variable 
(1) 


Annealed: 


N’ 


Note: indicates static collapse loads. 


were prepared for through Case-4 annealed specimens. The results are 


The dimensions such thicknesses, diameters, and lengths chord and 
brace are not always equivalent. Among these dimensions, there small 
variation diameter length brace thickness, they have nothing 
with the collapse load the specimen. With regard chord thickness, however, 
relationships between the collapse load and the chord thickness have been 
provided the form experimental equations. the present test, chord 
thicknesses were measured eight points each specimen, and their influence 
the collapse load evaluated. 

The variation stresses the neighborhood the Case-3 
specimen under static tensile load was analyzed finite element method where 
there was 10% variation the chord thickness. 

One eighth the specimen was divided into isoparametric shell finite 
elements with eight nodes. The circumferential stresses chord and the axial 
stresses brace are shown Fig. which the stresses corresponding 
10% decrease chord thickness are also indicated. Percentages increase 
stress caused 10% decrease chord thickness are represented Fig. 
The average increase stress the joint surrounded dotted line 


Data 
Normal: S, in kilonewtons 343° 392 294 304 275 265 
S”, in kilonewtons | 342* 277 301 295 272 264 
s"/ 1.024 0.943 0.895 0.892 0.814 0.787 
N 1 1 45 1 318 280 
N’ 1 87 141 176 295 306 
Ft S, in kilonewtons 334° 314 294 291 275 
S™, in kilonewtons | 340* 321 299 298 282 
i } | 0.990 0.929 0.868 0.865 0.814 
1 176 234 1 503 
q || 18 113 203 214 356 
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Fig. 10.2%. According this numerical result, concluded that 10% 
decrease chord thickness equivalent 10% increase stress, and this 
relationship would expected hold for moderate changes chord thickness 
and elastic behavior. 

the other hand, two kinds experimental equation have been reported 
for and K-type tubular joint specimens (4,7), which the static collapse 
load was proportional the chord thickness the 1/1.6 power 1/2 power. 
They are interpreted decrease chord thickness equivalent 16.5% 
21% increase static collapse Based both the numerical result 


FIG. 3.—Stress Distributions Chord and Brace: Stress Distribution 


and the experimental equations, the writers use the relation that the static collapse 
load proportional the chord thickness the 1/1.6 power, and assume 


the cyclic collapse load may proportional the chord thickness the same 
power. 


Strong scatter was observed the data number cycles failure 
The variation data was large that relationships among the data were difficult 
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FIG. 4.—Increase Stress Reducing Chord Thickness 10% 


clarify, and the writers introduce statistical approach based the Bayesean 
theorem. 

amount scatter even the case the simplest specimen, such smooth 
rod, because the fatigue phenomenon kind random process based 
the microcrack mechanism. However, the scatter the X-type welded 
tubular joint specimen described herein was large that cannot explained 
adequately ordinary fatigue mechanism. Number cycles failure 
seems strongly influenced welding defects, and the scatter 
increased when the scatter welding defects large. 

The scatter mainly result welding defects such lack weld 
penetration, condition weld surface smoothing, insufficient weld effective 
length, and degeneration material properties due welding. But mentioned 
the foregoing, since the weld line long and its geometrical pattern 
complicated X-type tubular joint specimen, difficult measure all 
factors (such weld penetration) along all weld lines, and impossible 
evaluate the influence various factors the cyclic collapse load. For 
example, close-up photographs the section including weld seam the 
neighborhood spot’’ the specimen are shown Fig. They show 
that weld penetration has been perfectly done, and the differences between 
them are difficult evaluate quantity. 


the only possible and practical resort, the following assumptions were 
made: 


Scatter welding defects can substituted scatter static collapse 


496 
¢ 
| 
A | brace 
1 
ele 
lo chord 
| 


WELDED TUBULAR JOINTS 497 


load, because static collapse load strongly influenced the welding defects. 

Scatter mainly result the difference static collapse load 
each specimen; that is, numbers cycles failure some specimens 
are extraordinarily small when the welding defects are very large (when the 
static collapse load very small), and vice versa. 

manufacturing error, errors cyclic load control, etc. 

there scatter static collapse load, that is, all specimens have 
the same static collapse load, the scatter dependent only the magnitude 
the repeated load. 


can seen Table the static collapse loads Case-3 specimens 
are scattered over wide range between 250 and 326 kN. Now 
let imagine that the fatigue test the repeated load 255 performed 
for these specimens. the case the weakest specimen whose collapse 
load 250 kN, the specimen bound collapse the first cycle 11), 
even before reaches the maximum repeated load 255 kN. However, 
the case the strongest specimen, whose collapse load 326 kN, the specimen 
may not collapse after hundreds load repetitions. Thus the scatter 
strongly and directly influenced the scatter static collapse load. 

order avoid this confusion, the modified value defined for 
the special case which the static collapse load not variable, that is, where 
there scatter defects. Let the static collapse load always 


FIG. Penetration near Hot Spot 
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equal certain value—an average (S,) used herein. the Case-3 
specimen, 290 shown Table The most probable value 

Assumptions Statistical Treatments.—The following three assumptions are 
introduced: 


The S-N curve straight line with negative gradient regular scales. 
high-cycle fatigue tests, almost linear relation generally obtained when 
the relationships and are plotted semilog log-log scales; but 
low-cycle fatigue test whose number cycles failure less than 10°, 
the situation not the same. According our test data including the static 
collapse load, approximate linear relation has been obtained when the data 
are plotted regular scales shown Fig. 2(a). The S-N curve defined 


which parameter determined from the data based the least square 
estimation. Parameter estimated herein —0.0175. 

becomes straight line with negative gradient. the high stress range, the 
plastic deformation will occur some parts the specimen, the stresses 
those parts are equalized, and finally the variation will reduce. 
consider the limited situation where the collapse loads are always equal 
standard deviation equivalent S,. The linear relation between 


which parameter determined from the data based the most likelihood 

The probability density corresponding arbitrary value repeated 
load assumed normal distribution generally done considering 
the probabilistic distribution (for instance, References and 5). The 
probability density collapse load also assumed normal distribution. 


Basic Concept Modification N.—The procedure modifying 
approach arranging random fatigue test data based the Bayesean 
theorem. The concept the approach explained the following referring 
Fig. 

The static collapse load the specimen, which used the fatigue 
test, unknown. Let consider that fatigue test performed for specimen 
whose collapse load unknown, under the repeated load magnitude 
and that the specimen collapses after cycles. The collapse load the 
specimen may S,, for example. These two situations are shown 
Fig. which the relation between and linear mentioned 
very small, because the situation located the tail part probability 
density for which called the other hand, the possibility 
other words, the collapse load the specimen S,, the specimen will 


as 
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FIG. 6.—Statistical Model 


seldom collapse the N*-th load cycle; whereas the collapse load S,, 
the specimen will possibly collapse the N*-th cycle. Both possibilities are 

Furthermore, the static collapse load follows its own probability density, 
which called Fig. 6(a), the possibility S,, which shown 

The conditional probability expressed, based the Bayesean 
theorem, 


P(S*, 
P(S,|S*, N*) = 


j 
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Since P(S*,N*|S,) equivalent and P(S,) equivalent 
F(S,)dS, and on, Eq. represented 


Here come the next stage, and the number cycles failure 
modified according the magnitude static collapse load. Let imagine 
that the specimen, whose collapse load S,, collapses the N*-th load cycle 
magnitude S*. have fixed the collapse load necessary 
calculate the number cycles failure corresponding under the 
condition that the collapse load This can realized equating the 
The value therefore obtained satisfy 


P(S,|S*,N*) 


(N) 


which the probability density g(N) normal distribution (Assumption 3), 
and can easily expressed 


Finally, the modified value N*, which called defined 
i 


which weighted mean the modified values according their 
probabilities occurrence. Eq. represented, using Eqs. and 


Test cycles failure Case-3 specimens are repre- 
sented Table between and nondimensional repeated loads 
are plotted the S-N diagram shown Fig. The condition 


the data before and after modification easily compared referring 
Fig. and (5). 


Size Errect ano Stress 


The influence relative size the specimen (size effect) will considered 
first, and after that the influence residual stress due welding, well 
the relationship between the two effects. The influence residual stress 
here called the annealing effect. 


WELDED TUBULAR JOINTS 501 


order evaluate the size effect, the static collapse loads and the numbers 
cycles failure under various repeated loads are examined using four 
kinds specimens different sizes, described the previous section and 
shown Table The size effect will evaluated quantitatively comparing 
the S-N curve Case-3 specimens with S-N curves Case-2, and 
Case-4 specimens. 

order evaluate the annealing effect and, the same time, the relationship 
between the size effect and the annealing effect, the collapse loads and the 
numbers cycles failure are examined using four kinds specimens 
different sizes, where half the specimens are annealed, shown Table 
The annealing effect will evaluated quantitatively comparing S-N curves 
annealed specimens with S-N curves nonannealed specimens, which are 
described specimens. should noted that there guarantee 
that the annealing effect really coming from the decrease residual stress 
successfully attained the annealing. could also consequence the 
change materials caused annealing. 

Test results concerning the size effect and the annealing effect are both 
represented Tables through which the scatter very large 
because the numbers specimens are few for both cases. The concept 
modification also introduced here order deduce S-N curve from 
random test data. 

Estimation S-N Curve.—First all, assumed that parameters and 
are always constant regardless size specimen and whether annealed 
not. According this assumption, the S-N curve estimated follows: 


Static collapse load first assumed, and then the S-N curve and the 
S-s, curve are fixed automatically. 
Repeated loads are modified due chord thickness, and numbers 
cycles failure are also modified assuming 
The sum the square distance between each data point and 
the S-N curve calculated. 


Both collapse load and S-N curve are determined minimize 
this sum. 


Size Effect.—The ratio the average collapse load Case-3 specimen 
290 KN) the estimated collapse load another case specimen 
Case-2, and Case-4 specimens are transformed the same level 
the Case-3 specimen multiplying The relation between 
and shown Fig. 2(a), and the relation between and shown 
Fig. 2(b). The S-N curve Fig. valid for every type specimen. 
Finally, let the reciprocal called r,. 

The ratio the cross-sectional chord area Case-2, and Case-4 
(A, 1913 called the size ratio r,, that is, A,/A,. The ratio 
r,, which called r,/r,, taken index size effect. 

weaker than the Case-3 specimen, vice versa. Parameters r,, r,, and r,/r, 
are shown Table and the relation between the size ratio and the index 
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shown Fig. which the relation obtained the least square 
estimation indicated thick line. Fig. ratio fatigue strength 
each specimen Case-3 specimen, for specimen; for 
specimen; and ratio cross-sectional area each specimen 


TABLE Effect and Annealing Effect 


Ratio 
Cross- 
Sectional 


Area Ratio Fatigue 


tional fa- resid- 
Area tigue ual 
square erty, prop- 
erty, 
(2) (8) (10) 


FIG. 7.—Size Effect 


Case-3 specimen. negative correlation between and r,/r, observed 
Fig. and concluded that the relative strength welded tubular 
joint specimen under low-cycle fatigue test shows downward trend when the 


q l 421 0.215 3.438 0.291 1.353 2.776 0.308 1.433 1.238 
i. 2 1218 0.591 1.593 0.628 1.063 1.481 0.577 0.976 1.076 
ra 3 1913 1.000 1.000 1.000 1.000 0.855 1.000 1.000 1.170 
iG a 3087 1.431 0.868 1.152 0.805 0.849 1.007 0.704 1.022 
1.6 
1.0 
> 
~ 
~ 
0.8 
7 
0.6 
0.0 0.5 1.0 
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size specimen increases. The reason for this could increase the 
welding defects, such lack penetration with increasing wall thickness, 
which would cause decrease strength. 

Annealing Effect.—The same process applied for the evaluation the 
annealing effect for the size effect. Instead the multiplier the ratio 
the collapse load the nonannealed Case-3 specimen the collapse load 
annealed Case-1, Case-2, and Case-4 specimens employed. This called 
the multiplier k’. The relations between and and between 
and are shown Figs. 2(a) and respectively. 

The ratio multiplier multiplier k’, which called k/k’, taken 
always constant for every case, there relationship between the size 
effect and the annealing effect. Parameters are shown Table 
and the relation between the size ratio and the index shown Fig. 

Two kinds straight lines can fitted Fig. They are: (1) dot-dash 
line which parallel the r,-axis; and (2) solid line which obtained 
the least square estimation (correlation coefficient 


FIG. 8.—Annealing Effect 


Following Line concluded that the strength welded tubular joint 
specimen increased more than 10% annealing. Following Line 
concluded that the size effect increased annealing. 


Low-cycle fatigue tests were performed for X-type welded tubular joint 
specimens. S-N curve was estimated using the modified values number 
cycles failure and after that, the size effect and the annealing effect 
were considered. They are summarized follows: 


Pull-pull tests, conducted for X-type welded tubular joint specimens, were 
restricted the low-cycle fatigue whose number cycles failure less 
than 600. 

The tests were performed sufficient number the same size specimens 
the case Case-3 specimens, convenient for the statistical treatments. 


The static collapse loads and the cyclic collapse loads were modified 
due chord thickness. 


1-4 
1-2 
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The most likely value was estimated for the restricted state, which 
the static collapse load every specimen equivalent. 

S-N curve was estimated using the modified value and 
(S” and where the S-N curve was assumed straight line regular 
scales. 

specimens different sizes were prepared order examine 
the size effect. The effect has been evaluated quantitatively based the index 
size effect (r,/r,). 

Annealed and nonannealed specimens were prepared order examine 
the annealing effect. The effect was evaluated quantitatively based the index 
annealing effect The tests were also performed four kinds specimens 
different sizes order evaluate the relationship between the size effect 
and the annealing effect. 

The following conclusions have been obtained: 


The statistical approach developed herein can applied the fatigue 
test data which are strongly influenced the scatter welding defects. 

case low-cycle fatigue, 600) welded tubular joint specimens, 
the S-N curve can assumed straight line the regular scales. 

Concerning the size effect low-cycle fatigue, the relative strength 
welded tubular joint specimens shows downward trend with increase 
the size specimen. 

Concerning the annealing effect low-cycle fatigue, the relative fatigue 
strength seems increase result annealing; but confirmed conclusion 


can obtained the relationship between the size effect and the annealing 
effect. 
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The following symbols are used this paper: 


parameters; 

chord diameter and brace diameter; 

probability density for 

probability density for 

multipliers; 

chord length and brace length; 

number cycles failure; 

modified value 

size ratio; 

1/k; 

index size effect; 

index annealing effect; 

repeated load; 

repeated load, which modified due chord thickness; 
tensile collapse load; 

collapse load, which modified due chord thickness; 
standard deviation 

chord thickness and brace thickness; and 

chord thickness (average the data). 


a,b,c 
d,,d, = : 
= 
N* = 
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BUCKLING WEB PERFORATED BEAMS 


INTRODUCTION 


Rectangular openings are frequently formed the webs steel I-beams 
facilitate the passage services, and have the effect reducing the plastic 
and buckling strengths the members. When two such openings are placed 
closely together, another potential mode buckling introduced, i.e., that 
lateral buckling the web post between the openings (10,13). 

This paper proposes method for determining the critical moment which 
buckling web post occurs. The validity the analysis demonstrated 
comparative results experiments scale models web posts; the 
experimental method described. The significance web post buckling for 
practical beams illustrated analyzing three different I-beams containing 
web posts varying widths. 


Collapse Modes.—There are two basic modes collapse for perforated sections 
containing adjacent web openings. the first these failure associated 
with only one hole, while the other involves the interaction both holes and 
the entire perforated section deforms single unit. Failure the web post 
prerequisite for the latter type collapse, and can occur either failure 
the ends the web post (Fig. 1), lateral buckling the post. Tests 
and 8WF20 beams (2,11) reported failure the ends the post, 
while buckling the post was reported test 14WF30 beam (10). 

Ultimate Load Analysis.—Fig. represents perforated section formed 
two adjacent concentric web openings length and depth separated 
web post width addition the primary bending moment, the 
section subject vertical shear force over its total length and thus 

Lect., Dept. Civ. Engrg., Univ. Natal, King George Avenue, Durban, 
Natal 4001, South Africa. 


Note.—-Discussion open until August 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication October 18, 
1979. This paper part the Journal the Structural Division, Proceedings the 
American Society Civil Engineers, Vol. 107, No. ST3, March, 1981. ISSN 
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has resist Vierendeel moment QL. For analytical purposes the section 
may represented the equivalent line structure Fig. where and 
are respectively the moments inertia the web post and the tee sections 
above and below the openings. The height the line structure taken 
the clear distance between beam flanges; this approximation was used 
Redwood (10), and its accuracy had been confirmed finite element stress 
analysis (2). 


(a) Dimensions (b) Line structure 


FIG. 2.—Perforated Section 


When collapse the perforated section involves failure the ends the 
web post, the resulting collapse mechanism contains six plastic hinges. This 
shown Fig. where only the upper half the mechanism shown—a 
plastic hinge occurs each end the tee section, and the web post opposite 
the tops the openings. 

The horizontal shear force the web post failure evaluated 
equating the moment F,(d/2) the plastic moment resistance the post; 


—— - 
FIG. 1.—Failure Ends Web Post 


ST3 PERFORATED BEAMS 509 


the retraining moment exerted the post the line structure thus F,(H). 
The ultimate resisting moment the entire perforated section Vierendeel 
bending therefore given 


which M,,, Mpc, and the plastic moments resistance 
the tees and respectively (2). 

When collapse the perforated section involves buckling the web post, 
the stresses the ends the post depend its slenderness and can vary 
from the fully elastic the fully plastic condition. the latter case the restraining 
moment the post identical that the post Fig. and the ultimate 
resisting moment the perforated section given Eq. When, however, 
these stresses are not fully plastic, the post will develop smaller restraining 


moment, and the ultimate resisting moment the perforated section will 
given 


which the horizontal shear force the web post buckling. The 
quantity represents the buckling moment the post, and general 


L/2 


Plastic 
hinge 


3.—Plastic Mechanism 


less than the restraining moment developed the fully plastic post. 
Therefore, determine the ultimate moment resistance perforated section 
such that Fig. means had developed computing the critical 
moment which the web post would buckle. 


The collapse mode steel beam subject the force system Fig. 
depends the slenderness the member. For short beam the ultimate 
flexural strength determined the full plastic moment the section, and 
failure due yielding the beam ends. For long beam the flexural strength 
determined the critical moment which the beam buckles elastically, 
and this depends the elastic stiffness the member. For intermediate 
span between these two extremes, failure occurs buckling after the beam 
has partially yielded. This yielding reduces the stiffness the affected regions, 
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and the beam’s buckling strength reduced accordingly. 

The analysis elastic beam buckling presents difficulty. The solutions 
number particular cases have been found solving the governing 
differential equations (15). Recent research using the finite element technique 
has provided more general solution the problem, and now possible 
(given suitable computer program) analyze elastic beam buckling for any 
arrangement loads, supports, and restraints 

contrast, the inelastic lateral buckling beams has not been extensively 
investigated. The majority early investigations had concentrated simply 
supported beams under uniform moment, and were directed towards the develop- 
ment theoretical model which took proper account the effects partial 
yielding the various factors controlling lateral stability. The problem the 
inelastic buckling single-span beams under moment gradient has been the 
subject recent studies (4,5,6), and these formed the basis the analysis 
Nethercot and Trahair (8) single-span beams under various load conditions. 


Full 
plasticity 


Inelastic 


buckling 


M=Fd/2 


FIG. 4.—Force System Web Post FIG. 5.—Inelastic Beam Buckling under 


Moment Gradient [Due Nethercot and 
Trahair 


The graph which they published for the analysis beams under moment gradient 
particular relevance the buckling web posts, and reproduced 
part Fig. 

Values the dimensionless inelastic critical moments were plotted 
against modified slenderness which inelastic buckling 
moment; plastic moment; and larger end moment elastic buckling. 
the original paper, curves were produced for the values 1.0, 0.5, —0.5, 
and —1.0. Only two these are shown Fig. i.e., 1.0, which corresponds 
the web post loading, and which the case uniform moment 
and given for purposes comparison. The difference the two curves 
illustrates the importance the moment gradient the inelastic range. The 
critical moments beams under uniform moment gradient —1) are reduced 
substantially below the full plastic moment, whereas those for beams under 


high moment gradient are reduced only slightly the point where 
full plasticity occurs. 


0.6 
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Evaluation M,, M,.—By evaluating and M,, Fig. can used 
determine the buckling moment single-span beam. The plastic moment 
given standard plastic theory, and the elastic buckling moment 
may found using equation due Clark and Hill (1), which applicable 
beams symmetrical section loaded either their ends, transverse 
loads applied the neutral axis. Thus 


rigidity; moment inertia about minor axis; torsion constant; 


coefficient; effective buckling length; and coefficient. 


For rectangular section such web post, and Eq. simplifies 


br 1/2 
(EI, GJ) 

Salvadori (12) has established the values the constants and for the 
type loading acting web post (Fig. 4). The value depends 
the degree which the beam supports restrain rotation about vertical axis 
through the beam ends, and varies from 0.5 for end providing full restraint 
1.0 for end permitting free rotation. usually assumed for design 
purposes that the web I-beam pinned the flange, i.e., This 
somewhat conservative assumption, since the flange restrains the web 
some extent. The present case complicated the fact that the degree 
restraint also depends the width the web post and the length and 
height the web openings; view the presence the latter, the usual 
value 1.0 will adopted for For K-value unity, varies between 
2.56 and 2.74, with the former value applying sections having negligible warping 
resistance; the value 2.56 will therefore used. 

Thus for web post width and thickness the elastic buckling moment 
can expressed 


2.567 


1.340 

Effective Buckling Length Web Post.—The calculation the elastic buckling 
moment web post complicated the fact that the post laterally 
restrained beam flanges situated clear distance apart [Fig. 6(a)], while 
the critical buckling stresses occur sections distance apart [Fig. 6(c)]. 
For the extreme case equal the critical moment given directly 
which the effective buckling length; for smaller values 
however, elastic buckling will occur larger moment. The problem may 


overcome assuming that this increase the elastic buckling moment 
effectively due reduction the buckling length the post. This 
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(a) Web post (b)Moments Stresses 


FIG. 6.—Web Post Moments and Bending Stresses 


reduced length will intermediate between and and may found 
comparing the buckling typical web post (Fig. with that post 
which Eq. applies, i.e., post uniform section having length and loaded 
Fig. 

The effective buckling length the uniform post and the critical stresses 
elastic buckling occur sections distance apart, thus 


and substituting from Eq. the horizontal shear force the post 
elastic buckling obtained. 


2.680 


P 
For the web post Fig. 6(a), the buckling length the post giving 
1.340 
but the critical stresses buckling occur sections distance apart [Fig. 
thus 


2.680 


therefore 


can inferred from Eq. that web post such that Fig. 6(a) 
may treated beam constant width and effective buckling length 
The elastic buckling moment then given 
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which true for post having zero depth and thus infinite width. 
was necessary however verify that Eq. was valid for intermediate values 
and this was done testing four scale models having range web 
post heights. 


INVESTIGATION 


The models Fig. were manufactured from single piece steel sheet 
and were proportioned give wide variation the height the web 
post; the choice 0.787 in. (20 mm) wide post ensured that elastic buckling 
would occur each case. Each model was subjected shearing forces 


1.97 high flange Triangular stiffener 
each side (except Model 


Plate thickness 0.046 
Note in. 25.4mm 


FIG. 7.—Model Dimensions 


12.60" long long 
free timber block fixed timber block 


Model clamped 
blocks 


Plate screwed 0.47" diam. steel rod 
block (far side also) 


(Clamps and clamping plates are omitted for clarity.) 
Note in.= 25.4mm 


FIG. 8.—Section through Test Arrangement 
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the plane the flanges reproduce the force system Fig. 6(a). 
The method applying the forces best described reference Figs. 
and 10. 

Each flange the model was placed contact with timber block. 0.236 
in. mm) thick mild steel plate was placed against the flange, and plate and 
block were clamped together (Fig. 9). One block was clamped the test table, 


FIG. View Test Arrangement 


FIG. 10.—View Connecting Mechanism 


- 
4 
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and will referred the fixed block. the underside the other (the 
free block) was screwed 0.118 in. mm) thick precision-ground steel plate, 
with the screw heads recessed into the plate. This rested similar plate 
placed the test table, and the contact surfaces were sprayed with dry 
lubricant just prior testing minimize frictional effects. steel mass 
was placed top the free block prevent uplift during loading (Fig. 8). 
With just this arrangement, the application horizontal force the free 
block would have caused rotate about the distant end the web post. 
was therefore necessary connect the two blocks such way that they 
remained parallel one another all times. longitudinal hole was drilled 
through each block, and 0.472 in. (12 mm) diam mild steel rod was inserted. 
The projecting ends the rod were threaded, and was secured position 
means square plate washers tightened against the ends the block. 
The extreme ends the rod were tongued, and were drilled receive steel 
pin. The elements linking the blocks comprised two mild steel bars 0.394 
in. (10 mm) square section; the ends the bars were grooved and drilled, 
and were pinned the tongued ends the rods form clovis joints. The 


TABLE 1.—Buckling Loads Model Web Posts 


inches inches inches Test Theory 


entire mechanism shown Figs. and 10; the latter figure the steel 
mass, clamps, and clamping plates have been removed for clarity. 

force was applied the free block through steel cable attached 
small welded steel A-frame; this frame was bolted the block means 
the threaded bar (Fig. 10), and was proportioned that the line action 
the force passed through the center the flange. The other end the 
cable was led over pulley clamped the end the test table, and was 
attached steel hanger which supported the applied loads. 

preliminary test was conducted establish the value the force required 
initiate sliding the free block; with the clamps and steel mass place, 
the magnitude this force was 6.75 (30 N). The models were then clamped 
turn the blocks, and load increments were applied until the web posts 
buckled; the view Model Fig. typical the deformed shapes 
obtained. Model was the only case where the web post received assistance 
from adjacent web material, and examination the deformed model revealed 
that despite the presence the clamping plates, the flanges had buckled locally 
their junction with the post. This had the effect reducing the buckling 
moment the web post. 


Model Theory 
number Test 
(1) (7) 
1.35 
4.57 5.19 126 1.09 
3.27 4.39 179 1.07 
1.97 3.41 287 
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The experimental shear forces causing buckling are given Table where 
they are compared with the theoretical values obtained from Eq. 11; the latter 
values were based elastic moduli 29.4 10° psi (203,000 
and 11.3 10° psi (78,070 N/mm’). 

The effect local buckling the flanges Model clearly shown 
the low test load, and this result should ignored. The use Eq. overestimated 
the elastic buckling moments the remaining models average 9%; 
Eq. could therefore regarded good approximation for the effective 
buckling length web post. 


using Eq. conjunction with Fig. was possible determine 
the buckling moment web post, and the state stress the critical sections. 
This approach was used investigate the significance web post buckling 
for practical beams. Three beam sizes were considered: (1) 8WF20 beam, 
which had the stockiest web the wide-flange beams available South Africa; 
(2) 16WF26 beam, which had the slenderest web any wide-flange beam; 
and (3) welded plate girder having in. (1,067 mm) deep 3/8 in. (9.5 
mm) thick web plate. The yield stress each web was taken 36,250 psi 
(250 and the elastic moduli and were taken 29.4 10° psi 
(203,000 and 11.3 10° psi (78,070 respectively. 

Each beam was considered perforated with pair concentric holes 


TABLE 2.—Values for Some Typical Beams 


Depth opening, Web post width, 
(2) (3) 
5.43 


CARN 


Plate girder 


Note: in. 25.40 mm. 


Beam 
(1) 
8WF20 
0.276 
0.339 
0.391 
0.479 
0.553 
10.43 0.387 
0.547 
0.670 
0.773 
0.865 
0.947 
1.023 
1.094 
28.00 0.602 
0.851 
0.951 
1.042 
1.203 
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near-optimum depth, i.e., two-thirds the overall beam depth the case 
the rolled sections, and two-thirds the web depth the case the plate 
girder. The quantity (M, was calculated for variety web post widths 
(Table 2), and the curve Fig. was used determine whether buckling 
would elastic inelastic, and whether the latter case would occur 
before the ends the post yielded completely. 

Buckling would not occur before the ends the post yielded fully, and was 
fact unlikely occur all. 

16WF26 beam.—For web post width in. (305 mm), Fig. predicted 
that buckling would occur just the ends the post yielded fully. For narrower 
widths the ends the post would fully yielded before buckling occurred, 
and buckling was unlikely occur all for the smaller widths. appeared 
that in. (356 mm) wide post would come very near buckling elastically, 
but when this width was compared with the 15.64 in. (397 mm) beam depth, 
was clear that the values for this and larger widths were 
academic interest only, since with strong post failure would confined 
single opening. 

Plate girder.—For widths in. (254 mm) and less, the ends the post 
would fully yielded before buckling occurred. For width in. (305 
mm) more, the post would buckle elastically. These values illustrate the 
narrowness the range between elastic and plastic buckling. 

These results show that for rolled sections, web post buckling will 
accompanied complete yielding the ends the post. The ultimate flexural 


strength the perforated section will therefore identical that given 
plastic analysis, and unnecessary perform buckling analysis the 
web post. plate girders, however, the slenderness the web may such 
that the post buckles while its ends are elastic only partially plastic. This 
will produce lower ultimate flexural strength than that given plastic analysis, 
and buckling analysis the web post mandatory. 


Summary 


The provision adjacent openings the web steel I-beam may cause 
buckling failure the post between the openings. this occurs while the 
ends the post are elastic only partially plastic, the ultimate flexural strength 
the perforated section will lower than that predicted plastic analysis, 
and buckling analysis the web post must performed. 

method has been presented for determining the buckling moment 
web post. expression was derived for the effective buckling length the 
post, and its validity was confirmed means tests scale models. The 
elastic buckling moment the post could then found from elastic buckling 
theory, and was used obtain the actual buckling moment from inelastic 
buckling analysis due Nethercot and Trahair (8). 

Analysis practical beams showed that web post buckling would not reduce 
the ultimate flexural strength perforated rolled sections; this was not the 


case with plate girders however, and with these structures buckling analysis 
the web post was mandatory. 
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The following symbols are used this paper: 


length web opening; 
width web post; 


4 
6 
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warping constant; 

depth web opening; 

modulus elasticity; 

horizontal shear force web post buckling failure; 
horizontal shear force web post plastic failure; 
modulus rigidity; 

effective height web post; 

moment inertia web post about major axis; 
moment inertia tee sections above and below openings; 
moment inertia web post about minor axis; 
torsion constant; 

coefficient; 

overall length perforated section; 

effective buckling length; 

inelastic buckling moment web post; 

larger end moment elastic buckling web post; 
plastic moment resistance web post section; 
plastic moment resistance tee sections; 

ultimate resisting moment perforated section Vierendeel 
bending; 

vertical shear force acting perforated section; 
thickness web post; and 

coefficient. 
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Basic ASPECTS BUCKLING 
COOLING-TOWER SHELLS 


Ali Yavuz Basar,’ and Ihsan 


INTRODUCTION 


The main loadings the hyperbolic natural-draft cooling towers are dead 
load and wind load. The biaxial stress state induced the axisymmetric shell 
through these loads nonaxisymmetric, changing along the height and circum- 
ference. Due this fact, the buckling problem becomes very complicated, 
even elastic range. 

There large number works the stability cooling-tower shells, 
both theoretical and experimental (4,5,6,7,9,10,11,14). The numerical results are 
very informative for axisymmetric cases (4,11). Under wind loading, however, 
the results obtained not seem realistic, due the nonperiodic buckling 
configuration this case According the observation made buckling 
tests under load control, the buckling initiated locally critical stress 
state reached (9). For this reason, Ref. critical stress state concept 
was proposed check the buckling safety cooling-tower shells. 

The objectives the present work are examine the effects boundary 
conditions, shape the meridian, and wall thickness. The parametric studies 
carried out yield numerical results which define some applicability limits for 
the buckling stress state concept for different cooling tower geometries. 


Linear 


detailed description the linear bifurcation theory used the present 
work given Refs. The equations equilibrium are derived considering 
the buckled state, which supposed infinitesimally neighboring. The 
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terms present these equations are considered consist two parts, one 
due the prebuckled state and another due buckled configuration. After 
eliminating the part due the prebuckled state, the following equations are 
obtained: 
+ + + 
+ 


aB a ' 


homogeneous part terms the terms the 
the bending membrane bending 
theory theory theory 


+ + 

The coefficients, L®, and the second strain tensor, allow the transition 
from the metric the deformed surface the metric the undeformed shell 
> a +, at + a + a 
+ + + + 

+ + + + + + 
The constitutive equations 


+ 1 


together with the symmetry condition 
+ 


and the Love-Kirchhoff-hypothesis 


complete the system the differential equations necessary determine the 
critical load parameters. Eqs. 1-10, the nonsymmetrical stress resultant 
tensor; stands for the symmetrical tensor the stress couples; and 
the shear force tensor. Superscripts and denote the prebuckled buckled 
state, respectively. The variable the elasticity tensor, and being 
the modulus elasticity and wall thickness the shell. Eqs. and 


the last terms are neglected, then the equations the so-called classical 
buckling theory are obtained. 


For general shells revolution, according Fig. the differential Eqs. 
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1-10 can written explicitly using the rules tensor calculus. the variables 
are given column vector 

The partial differential equations the problem can written matrix 
notation as: 


+ 
homogeneous part stability terms 
the bending 
theory 


stand for the partial derivatives with respect the cylindrical 
coordinates and respectively, according Fig. 

For axisymmetric prebuckling stress state with 
nr? =g'=m"=0 
the partial differential equation system can reduced system ordinary 
differential equations, for the symmetric and antisymmetric variables, the 


xis 


x! 


FIG. 1.—Curvilinear Coordinates for Shells Revolution 
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following statements are introduced: 


The variables and are now only functions order eliminate 


the derivatives higher order new variables are introduced: 


+ 
and new vector the unknowns obtained: 
+ toy toot) to + + + 
T i 12 22 
(n,, n,, n,, n, GWnIn™, m,, m, W., 
WinWrn Sn S2n Vin Von Van) 


Written matrix notation, the final system ordinary differential equations 

+ + + 
AY’ 


The variables and are square matrices the size which depend 
the number waves assumed along the circumference. These matrices 
are given Ref. The variable the critical load parameter. 

For integration, the differential equations are transformed into finite difference 
equations using the Hermitian (multilocal) method. The shell divided into 
equal intervals and the corresponding finite difference equations are 
written applying the Simpson’s rule for the intermediate intervals and the Bessel’s 
formula for the first and last intervals follows: 


ith 


(1) (2) 


Applying the summation operator Eq. 19, and considering Eqs. and 
21, equations the internal nodes (2) are obtained 


3 + + 3 + 


The equations the first and last nodes (1) and (k) are derived similar 
way using Eqs. and 23. Afterwards, the critical load factor, calculated 
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solving the eigenvalue problem the form 
K,X + Aner 0 


The boundary conditions, four each boundary, are considered, easily substitut- 
ing them the corresponding equations, explained Ref. 


The experimental results considered this work are those Ref. 
the tests described detail Ref. the buckling procedure was always initiated 
locally, the experimental values for the symmetrical models, and for 
being the same range. other words, the tests showed that the reduction 
the height above the throat from 600 (23.65 in.) 220 (8.65 in.) 


FIG. 2.—Power Versus Ratio 


doesn’t affect the buckling stresses and the buckling behavior. This fact can 
explained through the buckling initiating effect the local imperfections, 
which was the same both tests series. addition, this initially surprising 
result agreement with the buckling behavior cylindrical shells medium 
length subjected axial compression. The buckling stress independent 
the shell length, both tests and according the theory, and does not change 
clamped hinged boundary conditions are present. 

The buckling loads obtained from the tests, given Fig. Ref. 
are extrapolated for wall thickness (0.0785 in.). For that purpose, 
accordance with the theoretical results obtained and given Fig. the 
present work, the buckling stresses measured are multiplied for 
each model. The results given Fig. are obtained this way, and are identical 
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The first letter, stands for the type the model symmetrical 
cooling-tower shaped. The second letter denotes the boundary conditions 
considered. The letter stands for clamped ends both boundaries. case 
the lower boundary the model assumed clamped. The upper 
boundary also clamped; however displacement the vertical direction 
allowed, indicated Fig. the tests, the upper head plate the 
testing machine was resting free the model, vertical displacement rotation 
being possible. However, different from case the upper edge the model 


KN/mr? (psi) 


| 
Model Number 
Boundary conditions 
Shape 4 
Experimenta! results 


o SM! 
x SM2 


| 


| 
| 
| 
| 
| 


FIG. 3.—Test Results and Theoretical Diagrams 0.0785 in. mm); 
500,000 psi (3,450 and 0.38] 


also stuck the head plate, and only displacements and rotations the relatively 
rigid head plate plane were possible. Therefore, the real boundary conditions 
the tests were between and 

comparison the interaction diagrams with the test results provides the 
following information. The calculation results different interaction diagrams 
for the symmetrical and cooling-tower models, fact contradiction the 
tests. The interaction diagrams and can nearly accepted the upper 
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and lower boundaries, respectively. For stress states with nearly zero stress 
the circumferential direction, the test results match with the values calculated 
for symmetrical models having the boundary conditions quite well. For other 
stress states the interaction diagrams fit better. 

The discrepancies Fig. can explained mainly means the assumption 
and approximation the buckling analysis used. For example, the initiation 
the buckling procedure the tests locally, whereas the theory assumes 
instantaneous periodic buckling pattern. Furthermore, the parametric 
investigations the following section shows, the buckling loads calculated are 
substantially influenced the boundary conditions. 


The buckling loads due lateral pressure and axial compression were obtained 
for the boundary conditions and according Fig. The calculated 


450mm 


FIG. 4.—Boundary Conditions Considered Parametric Investigation 


shells have symmetrical shape with respect the throat, the meridian curve 
being circular arch. The height the shell, and the radius the base, 
R,, are 600 (23.6 in.) and 225 (8.85 in.), respectively. The radius 
the throat, varied steps 0.025, the range 0.60 
0.85. The wall thickness, 1.6 (0.063 in.), and material 
constants, 3650 (530.000 psi), and 0.35 are considered. 
applying the finite difference method, the meridian divided into equidistant 
intervals. 


Fig. the buckling loads the symmetrical shells with boundary conditions 


| 
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and under lateral radial pressure are given for different ratios 
with the corresponding number buckling waves, the circumferen- 
tial direction. The boundary conditions, are nearly the same those present 
natural-draught cooling towers, whereas the type supporting, corresponds 
the so-called classical boundary conditions. Boundary conditions can 
considered representative for cooling tower ground stressed 
the yielding and deforming the direction the supporting columns without 
any resistance, hypothetical lower limit. 

The curves for the buckling loads change regularly for system SB, SF, and 
SI, the buckling pressure decreasing with increasing ratio. Also, the 


R, /Re 


FIG. 5.—Buckling Pressure Versus Ratio 


number the buckling waves changes and slightly. contrast 
these systems, has buckling loads and wave numbers alternating strongly 
due slight variations This high geometric dependence can hardly 
accepted realistic. hyperboloidal shell having classical boundary conditions 
tends deform quasi-inextensional, i.e., eliminating the membrane energy. This 
behavior results low buckling loads, shown Fig. for the cases 
0.75 and 0.85, with and the other hand, for 0.65, 
the deformation mode with five periodic waves the circumferential direction 
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comprises considerable extensional part. deformation mode combination 
modes with different wave numbers, might comprise much lower extensional 
deformations, and might yield lower buckling load. However, the present 
analysis does not include solutions with buckling forms combination 
different wave numbers. For this reason, approximation, the envelope 
the buckling curves, which again regular curve, can considered 
realistic for shells having classical boundary conditions. The boundary conditions 
allow the range 0.60 0.85 overall inextensional deformations 
with For this reason, the buckling curve, very low and 
regular. 


FIG. 6.—Buckling Axial Load Versus Ratio 


The buckling loads under axial compression are given Fig. Again, for 
the boundary conditions and the curves are regular, having increasing 
tendency with increasing ratio. For the shell having classical boundary 


conditions, the same interpretation the case external lateral pressure, 
valid. 


The shape the meridian natural-draught cooling-tower shells mostly 
hyperbola. For this reason, the parameter studies were extended also 
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FIG. 8.—Buckling Axial Load for Shells with Circular and Hyperbolic Meridian 
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such shells with the boundary conditions which are representative for cooling 
towers. 

Fig. shows the critical lateral pressure function the ratio 
for the cases the circular meridian, SF, and the meridian hyperbola, 
SFH. The fact that the critical lateral pressure is, case SFH, almost independent 
the ratio, very perplexing. addition, the divergency the curves 
and SFH with decreasing ratio may unexpected, due only 
slight difference the shape cooling tower with circular hyperbolic 
meridian, the same for both cases. interpretation this initially 
contradictory behavior possible the buckling tests and buckling stresses, 
Ref. are considered. 

tests models with circular meridian, the buckling always initiated 
the throat (9). This due constant curvature along the meridian. the 
case hyperbola the curvature decreases with increasing distance from the 
throat. For this reason the buckling initiating point tends remove regions 
away from the throat, and the first buckling wave forms near the boundaries. 
This tendency was observed tests models with hyperbolic meridian 
(11). Due bigger radius curvature, the circumferential stress initiating 
the buckling will also larger than the stress the throat, under the same 
lateral loading. Fig. the ratio the buckling loads for SFH and 
approximately equal the ratio. That means the buckling initiating 
stress remains the same for shell revolution having negative Gaussian 
curvature, for given and R,, independent the shape the meridian. 
This explains why operate with critical stresses stress states instead 
critical loads (9,10). 

The buckling loads the case axial compression are given Fig. 
The curves have the same tendency both for system and SFH, the values 
for the case hyperbola being larger. Under axial loading, the biggest meridional 
stresses are acting the throat. For this reason, obvious that models 
with circular meridian buckling always the throat (9) can carry lower axial 
load than the hyperbolic models buckling regions away from the throat (11). 
the stresses are concerned, turns out, that both cases the buckling 
initiating point, approximately the same critical stress interaction relationship 
Ref. holds. This fact allows the extension the equations Ref. 
shells with hyperbolic meridian and regions other than the throat. 


addition shells with 1.6-mm (0.063-in.) wall thickness shells with 0.8-mm 
(0.0315-in.) and 0.4-mm (0.016-in.) thicknesses were also calculated. this way, 
range shell slenderness 1/565 1/85 for 0.60 1.00, 
could covered. Circular meridian and the boundary conditions System 
were considered. The calculations were carried out separately under lateral 
pressure and axial compression. 

The results this parameter study are given Fig. expressions for 
the critical loads, and the wall thickness the shell separated and 
considered the slenderness nondimensional form, with being 
the radius the throat. The and functions the other geometric 
parameters and material properties. Fig. shows, both critical lateral and 


MARCH 1981 ST3 


loads are proportional with nearly constant power the wall thickness 
the range 0.60 0.90 which interest for hyperboloidal 
cooling-tower shells. For cylindrical shells, the well-known factors 5/2 
and are obtained. 


The following conclusions can made: 


The agreement between theoretical results and experimental results for 
hyperboloidal shells better than the case cylindrical shells. For this reason, 
the results the bifurcation theory for cooling-tower shells under axisymmetric 
loading can accepted realistic and used for design after slight correction. 

upper cornice ring, reducing the displacements normal and circum- 
ferential directions the top the shell, proves advantageous, the 
comparison the buckling curves for systems and shows. 

For the boundary conditions which are general identical with those 
cooling-tower shell, the relationship between the buckling loads and 
regular, and important geometric parameter the cooling-tower 
buckling. 

the case settlement the supports the meridional direction, the 
buckling resistance the cooling tower shell drops considerably, the results 
the limit case boundary conditions demonstrate. 

Although different buckling loads are obtained for hyperbolic and toroidal 
cooling-tower shell having the same ratio, the difference vanishes 
the buckling stresses the points initiating the buckling process are considered. 
That important reason why the buckling design should based the 
stresses and not the loads. 

For cooling-tower shells, the buckling loads, meridional compression 
uniform lateral pressure, are proportional Therefore 
the parameter taking account the slenderness the shell the expressions 
Ref. for the buckling stresses. 
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Square matrices (20 20); 


b,b,,b, square matrices (14 14); 
tensors; 
geometry and material properties for axial 
load and external pressure, respectively; 
height the shell; 
i,,i,,i, unit vectors; 
K,,K, matrices; 
number modes the finite difference method; 


transformation coefficients; 


tensor the stress couples; 
power for axial compression and lateral pres- 
sure, respectively; 
shear force tensor; 


square matrix (20 20); 


te, 

Wes 

Y,Z 
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new variable; 

summation operators; 

square matrices (14 14); 

wall thickness shell; 

critical wall thickness shell; 

new variable; 

displacements the shell middle surface; 

solution vector; 

cartesian coordinates; 

vectors the unknowns; 

index 

first strain tensor the shell middle surface (extensions); 

curvilinear coordinates; 

critical load parameter; 

ratio; 

stresses circumferential and meridional directions, 
respectively; 

second strain tensor (change curvature); 

prebuckled and buckled states, respectively; 

partial derivatives with respect and respectively; 
and 

covariant differentiation with respect 
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SIMPLIFIED ANALYSIS EDGE STIFFENED 
CANTILEVER SLABS 


Baidar Bakht' 


INTRODUCTION 


Some components bridge and other structures can designed treating 
them cantilever slabs isolation without any appreciable loss accuracy. 
For the purpose analysis these components are usually assumed composed 
homogeneous isotropic material exhibiting linear elastic characteristics. When 
the cantilever slab subjected concentrated loads accurate assessment 
the force effects could often obtained either computer based methods 
(e.g. grillage analogy finite element methods) rigorous methods 
Jarmillo (5) and Silberstein (12)], that are too complex for every day design 
office applications. 

Simplified methods, suitable for design office use, are available for cantilever 
slabs linearly varying thickness but having edge beams (2,3). One possible 
way using these simplified methods for the analysis cantilever slabs with 
edge beams replace the edge beam slabs equivalent flexural rigidity 
and then analyze the resulting slab cantilever slab without edge beams. 
Such approach produces satisfactory results only the flexural rigidity 
the edge beam small compared that the slab. When the edge beam 
rigidity large, the equivalent slab approach does not provide results acceptable 
accuracy. 

Fig. shows the distribution longitudinal moments along the cantilever 
root due concentrated load, cantilever slab with edge beam, and 
slab without edge beam. The moments were obtained 
analyzing the two slabs the grillage analogy method which examined 
later the paper. can seen Fig. that the equivalent slab approach 
fails correctly reproduce the behavior the edge-stiffened slab. 
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box-girder bridges those free edges the cantilever slabs that are parallel 
the span the bridge are often stiffened beams. These beams are usually 
provided for functions other than the enhancement the behavior the slab, 
such providing curbs barrier walls for the road. Even though their presence 
may incidental, the edge beams enhance the transverse distribution 
concentrated loads cantilever slabs extent which makes their consideration 
the analysis quite prudent. 

simplified semigraphical method for determining longitudinal moments 
edge stiffened cantilever slabs linearly varying thickness given this 
paper. The paper also contains methods for determining maximum moments 
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70 


DISTANCE ALONG THE ROOT FROM LOAD 
FIG. Slab Method Analyzing Edge Stiffened Cantilever Slabs 


the edge beam due concentrated loads anywhere the slab. Other responses 
such deflections and beams cannot obtained the proposed methods. 
The methods given the paper are based results obtained analyses 
using the grillage analogy method. Therefore, the proposed methods are subject 
the same limitations and errors the grillage analogy method which briefly 
examined the following. 


Lightfoot (6), Yettram and Husain (13), and Sawko (8) are among many who 
have suggested that the behavior isotropic slab can predicted 
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equivalent assembly orthogonally connected beams, which can solved 
any one the (by now) many commercially available grillage programs. 
The writer (1) has shown that for slabs with zero Poisson’s ratio all the proposed 
idealizations converge common approach which the flexural rigidity, 
EI, beam the equivalent grillage equal the flexural rigidity 


2.0 


GRILLAGE ANALOGY 


SOLUTION) 


My, 


FIG. 3.—Plan Infinitely Wide Cantilever Slab 


the strip slab that represents. The torsional rigidity, GJ, the beam 
taken numerically equal its flexural rigidity. This concept idealization, 
which was proposed for slabs uniform thickness, could also extended 
slabs linearly varying thickness which case the idealized grillage consists 
beams varying rigidity one direction and those constant rigidity 
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the other. all points the torsional rigidities beams should equal 
their corresponding flexural rigidities. The grillage programs for analyzing 
the resulting grillage should capable dealing with beams varying properties, 
such the program described Ref. 11. The aforementioned program based 
techniques described Sawko and Wilcock (10). 

The validity the grillage analogy for cantilever slabs linearly varying 
thickness demonstrated Fig. which shows the distribution longitudinal 
moments for various load cases along the cantilever root very wide cantilever 
slab linearly varying thickness. The results obtained the grillage analogy 
method are compared with those given Homberg and Ropers (4). The latter 
results were obtained using classical plate theory assuming small deflections 
and shear deformations. Both analyses assumed the Poisson’s ratio 
equal zero. From the results plotted Fig. can seen that the grillage 
analogy gives results identical those predicted the plate theory. 

The assumption zero Poisson’s ratio necessary able idealize 
the slab orthogonally connected beams. Consideration finite value 
Poisson’s ratio necessitates the introduction disposal members the grillage 
(13). Neglect Poisson’s ratio results underestimation the longitudinal 
moments. The underestimation, however, usually small being less than about 
concrete slabs. 


The simplified solution proposed the writer and Holland (3) for cantilever 
slabs with linearly varying thickness but without edge beams given the 


form 
PA’ 


which the longitudinal moment; and constant given graphical 
form for various slab thickness ratios, distances the load, and reference 
points from the cantilever root. The notation for other variables shown 
Fig. The solution, which modelled previous solution Sawko and 
Mills (9), provides results comparable those obtained rigorous methods; 
satisfies the following essential conditions: 


The area under the longitudinal moment curve plotted for given value 


The introduction edge beam cantilever slab does not change any 
the foregoing conditions. Therefore seemed probable that the solution 
the writer and Holland (3) with different values could also applicable 
edge stiffened cantilever slabs infinite width. The main problem with 
such proposition remained that characterizing the edge beam parameters 
nondimensional form. literature search revealed rigorous solution 
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Reissmann and Cheng (7) for edge stiffened cantilever slabs constant thickness 
and infinite width. The series solution given terms two nondimensional 
parameters, k,, and k,, given 


which and GJ, the flexural and torsional rigidities the edge beam; 
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which the moment inertia the longitudinal section the slab 
about its middle surface. 

Reismann and Cheng (7) show that the behavior edge reinforced slab 
affected mainly and that changes have relatively little effect 
the manner which deflections and longitudinal moments are distributed. 
For example, they show that very wide cantilever slab uniform thickness 
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the maximum longitudinal moment the cantilever root due concentrated 
load the free edge increase only 0.34% due 400% increase k,. 
seemed reasonable delete from consideration and assume that the load 
distribution behavior very wide cantilever slab can characterized 
only k,. Later calculations, which are not reported this paper, justified this 
assumption. 

Fig. demonstrates that one the parameters characterizing the 
distribution longitudinal moments edge stiffened cantilever slabs; shows 
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the distribution cantilever slabs different dimensions but having 
the same value k,. can seen that the pattern distribution 
the same the two slabs. 

Several analyses using the grillage analogy method confirmed the validity 
Eq. for cantilever slabs linearly varying thickness and having edge 
beams. The values coefficient for given values and can obtained 
from the corresponding maximum values given the grillage analogy method, 
from the following equation: 


(M,), =-0 = 


Once the value known the cantilever slab can analyzed using 
Eq. selection the several comparisons between the grillage analogy 
results and the proposed solution given Fig. can seen that Eq. 
can predict the values close conformity with those given the 
grillage analogy method. 

Cantilever slabs having values and and having edge beams 
with equal 7.5, and were analyzed the grillage analogy method 
for various load positions. take advantage the symmetry the structure 
about the load only half the slab was analyzed. The edge the half structure 
that coincides with center line the full structure was prevented from rotating 
about the x-axis thus simulating the symmetry. Along the y-axis the length 
the plate was taken equal 4a. The grillage mesh employed was 
the same that used the writer and Holland (3). Values obtained 
from the foregoing analyses are plotted Fig. which conjunction with 
can now used analyze wide edge stiffened cantilever slabs provided 
that the concentrated load least distance away from the free edges 
parallel the x-axis. 


Most cantilever slab components box-girder bridges have unsupported 
longitudinal edges the vicinity which concentrated loads are indeed applied. 
order completely analyze the slab, there should method comparable 
simplicity that could applied the slab regions near the longitudinal free 
edges. The writer al. (2) have provided the following solution for 
the root semi-infinite cantilever slabs without edge beams, which one 
longitudinal free edge near the concentrated load but the other least 
distance from the load: 


which the same coefficient that for infinitely wide cantilever slabs, 
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examined earlier; coefficient which depends upon and the load 
position with respect longitudinal free edge and the transverse free edge; 
other notation shown Fig. The foregoing solution, which satisfies 
all the necessary boundary conditions, has been shown the writer, al. 
(2) provide solutions comparable those predicted computer based rigorous 
linear elastic methods. 

noted that Eqs. and were not rationally derived but chosen satisfy 
the necessary boundary conditions and match the rigorously obtained solutions 
most closely. Eq. the first term the parenthesis provides solution for 
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FIG. 10.—Edge Beam Moments 


the infinitely wide plate and the second term represents the increase the 
longitudinal moment that induced due the introduction longitudinal 
cut the infinitely wide slab. 

Several analyses, similar the earlier examined analyses for wide cantilever 
slabs, showed that Eqs. and which were developed for cantilever slabs 
without edge beams, were also applicable edge stiffened cantilever slabs 
with good degree accuracy. selection comparisons between the grillage 
analogy results and those obtained and given Fig. For 
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Thus from known values the cantilever root the longitudinal free 
edge and A’, the value can obtained from Eq. 

large number edge stiffened cantilever slabs, having various values 
and were analyzed for different load positions the vicinity one 
the longitudinal free edges. The length the slabs along the y-axis was 
kept sufficiently long keep one free longitudinal edge least distance 
away from the load. From the results the analyses and already developed 
values A’, values coefficient were calculated using Eq. the resulting 
values are plotted Fig. Using the relevant value from the 
graphs and the value from Fig. cantilever slab with concentrated 
loads the vicinity free longitudinal edge can manually analyzed 


Beams 


When concentrated load applied edge stiffened cantilever slab, 
the edge beam subjected both sagging moments (causing tension bottom) 
and hogging moments. addition being influenced values and 
the magnitude moments edge beam also governed the length 
the cantilever. was found that the nondimensional quantity (edge beam 


Pa) remained constant slabs having the same values and 
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but different lengths. Fig. 10(a) shows the distribution (edge beam moment Pa) 
two cantilever slabs infinite extent. Both slabs have 2.0 and 
7.5, but one slab has cantilever length units and the other 160 
units. 
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FIG. 12.—Values Coefficients and 


can seen that inspite the large difference the cantilever length 
the manner distribution edge beam moments very nearly the same 
the two cases. The foregoing exercise was repeated for cantilever slabs 
semi-infinite extent, Again was found, shown Fig. that the 
nondimensional quantity (edge beam moment/Pa) was the same the two 
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cases. was concluded that the maximum edge beam moments could obtained 
using the following nondimensional coefficients: 


coefficients for the maximum sagging and maximum hogging 
moments the edge beam infinitely wide cantilever slab, respectively. 

coefficients for the maximum sagging and maximum hogging 
moments the edge beam cantilever slab semi-infinite width, respectively. 


The maximum hogging sagging moment the edge beam equal the 
product and the relevant coefficient. Values the various coefficients 
were obtained from the already conducted analyses that were done get the 
values and B’; these values are plotted Figs. and 12. The maximum 
sagging moment the edge beam always occurs (for notation see 
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FIG. 13.—Values Y,/a 


Figs. and 7). The maximum hogging moment occurs distance, Y,, from 
the load along the y-axis. The value was found relatively independent 
series curves were developed from the results the grillage analogy analyses 
obtain fairly realistic estimates Y,; these curves are given Fig. 13. 
reminded that these curves are intended provide only fair estimates 
with accuracy about +10%, and not exact values Y,. 


edge stiffened cantilever slabs the longitudinal moments anywhere the 
slab can obtained Eq. conjunction with Fig. the concentrated 
load least distance away from both the longitudinal free edges. When 
the concentrated load the vicinity longitudinal free edge the longitudinal 
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moments along the cantilever root can obtained Eqs. and conjunction 
with Figs. and Locations and values maximum moments edge stiffened 
cantilever slabs infinite and semi-infinite width can obtained with the 
help Figs. 11, and 13. 


Simplified semigraphical method for manual application are presented for linear 
elastic analysis edge stiffened cantilever slabs. The methods, which are based 
previous solutions for cantilever slabs without edge beams, require very 


little effort and yet provide results comparable those obtained computer 
based rigorous methods. 
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The following symbols are used this paper: 


dimensionless constant defined Eq. and obtained from Fig. 

length cantilever slab measured perpendicular supported edge; 

dimensionless constant defined Eq. and obtained from Fig. 

distance concentrated load from supported edge measured along 
x-axis; 
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flexural rigidity slab per unit length; 

distance concentrated load from longitudinal free edge, measured 
along y-axis shown Fig. 

modulus elasticity; 

coefficient for maximum sagging and maximum hogging moments 
edge beam infinitely wide cantilever slab, respectively; 
coefficients for maximum sagging and maximum hogging moment 
edge beam cantilever slab semi-infinite width, respectively; 
shear modulus; 

moment inertia grillage beam; 

moment inertia edge beam; 

moment inertia longitudinal section slab about its middle 
surface; 

torsional inertia grillage beam; 

torsional inertia edge beam; 

torsional inertia longitudinal section slab; 

coefficient defined Eq. 

ratio flexural rigidities edge beam and cantilever slab; 

ratio torsional rigidities edge beam and cantilever slab; 
longitudinal moment per unit width; 

intensity concentrated load; 

slab thickness free edge; 

slab thickness supported edge; 

distance measured perpendicular and from supported edge; 
distance maximum hogging moment from concentrated load, 
measured along y-axis; 


distance measured along supported edge from concentrated load; 
and 


distance measured along supported edge from longitudinal free 
edge, shown Fig. 
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INTRODUCTION 


Diaphragms are common and economical structural systems for resisting lateral 
forces that develop from wind and earthquake. typical diaphragm consists 
several layers beam elements that are bonded together shear connectors. 
Loadings act the plane the diaphragm. The system highly indeterminate 
and its analysis complex (Fig. 1). 

The layers diaphragm are connected together shear connectors. The 
connector may nails, glue, bolts, screws, welded spots, 
shear clips. Regardless the shear connector used, the connector deforms 
when the diaphragm loaded and allows the layers slip. Thus partial 
composite action between the layers developed. 

The design diaphragms complex problem and, present, employs 
empirical methods which are based experience and testing. 

The first writer and Hiremath (8) presented theoretical method for the analysis 
diaphragms. The method develops set differential equations for the 
analysis multilayer diaphragm. The solution these equations was obtained 
using finite difference method. The study, however, did not present experimental 
verification results. 

Almost all the lumber diaphragms will have courses containing members 
different lengths make the full length the diaphragm. This causes 
discontinuities the courses. addition, diaphragms may have boundary 
members, with effects not included the theory that has been presented 
the first writer and Hiremath (8). 

The purpose this paper present experimental comparison with 
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the theoretical procedure Ref. This accomplished constructing 
lumber diaphragm without any discontinuities boundary members and subject- 
ing uniform loading. Axial strains the different layers and deflection 
the midspan are measured. 

second objective this study observe the effects discontinuities 
and boundary members diaphragm behavior. evident that the discontinui- 
ties the layers will permit increased deflection and strain, while the boundary 
members will reduce these effects. Magnitudes these effects are evaluated 
from this experimental study. 

The third objective evaluate the effect connecting joists the diaphragm 
with glue. This will develop resisting couples that will add the stiffness and 
strength the diaphragm. 


Review 


brief review the major previous development leading the study 
diaphragms presented here. The study diaphragms with layers can 


w ibs/ft 


FIG. 1.—Typical Diaphragm System 


classified broadly the area partial composite action incomplete interaction 
between structural components. 

The basis theoretical study the area incomplete interaction the 
pioneering work Newmark (13) 1951, who studied the incomplete interaction 
composite T-beam having concrete deck connected steel I-beam 
shear connectors. formulated second-order differential equation 
express the relationship between slip and the longitudinal force 
each element the T-beam. few years later, Clark (2) reported his study 
laminated beams where employed the moment-area method compute 
the deflection beam with two more laminations. Goodman (4) extensively 
studied the interlayer slip problem involving layered systems. 

The third writer and McGee (7), using the same type differential equation 
Newmark, analyzed wooden T-beams consisting wood joist and deck 
connected with elastomeric adhesives. They considered concentrated load 
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any point the beam and distributed load. For these loading conditions 
they obtained the closed form solution this system. Using the superposition 
method, solutions for other types load are possible. Similar solutions were 
also reported Vanderbilt (15). Expansion the foregoing approach 
three-layer system subjected uniformly distributed load was done Anderson 
(1). The first writer and Brito (9) developed closed form solution for two-layer 
system with gaps one layer. The interlayer slip this study applicable 
both adhesive bonds and mechanical connector links the compliance 
relationships can determined. The only difference that the case 
adhesives the width and thickness the glueline must accounted for while 
the mechanical connectors the spacing must considered. 

the number layers becomes greater than three, becomes very difficult 
use closed-form solution. Finite difference and finite element techniques 
may used. (12) applied finite difference technique solve the governing 
equations for composite beams having axis symmetry and arbitrary number 
layers fastened together with mechanical connectors. Foschi (3) developed 
computer program analyze plywood diaphragm using finite element 
techniques. His study accounted for the orthotropic plate action and nonlinear 
connection behavior. The analysis gives good estimates for the deformation 
diaphragms and ultimate loads related connection yielding. The first writer 
and Hiremath (8) extended the three-layer system solution developed Anderson 
(1) system with arbitrary number layers using finite difference 
solve the equations. 

Diaphragms for floor and roof construction have been tested the Forest 
Research Laboratory Oregon State University, the American Plywood Asso- 
ciation, and the United States Forest Products Laboratory develop data for 
the design buildings. Wood diaphragms are usually sheathed with either plywood 
lumber. Johnson (10) tested two diaphragms in. in. the year 
1968. One used nailed connections only and the other employed nails and 
elastomeric adhesive provide additional strength. found that the roof 
section with adhesive was several times stiff similar roof section without 
adhesive. The diaphragm with adhesive sustained three times the maximum 
load the nailed system. also made some openings the diaphragms 
and studied the effect performance. found deflection increased due 
the openings the glued diaphragms and the recovery the diaphragm after 
the release this load was retarded. Later Johnson tested four more diaphragms 
the same size (11). One these was 100% glued, the second one was 50% 
glued, the third one was 100% mechanically fastened, and the fourth one was 
50% mechanically fastened. The fasteners were shear clips. was found that 
maximum loads for roof sections with 100% shear clips adhesive were about 
1.6 times greater than maximum loads for corresponding roofs with 50% shear 
clips adhesive. Diaphragms with shear clips were about 15% stronger than 
those with the elastomeric adhesive used this study. All four roof sections 
were two three times strong simple nailed diaphragm with neither 
shear clips nor adhesive (this result would greatly affected frequency 
nailing and nail size). also found that shrinkage wood due change 
moisture content caused about the total length glue line 
cracked. This directs attention problem glue systems which merits some 
consideration. 
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brief review the theory presented the first writer and Hiremath (8) 
examined next. 

The partial composite action interface connected layers develops 
shear force. This shear force can represented (refer Fig. couple 
acting the centroids the two layers. The sum the moments acting 
the two layers given by: 


external moment section the beam; the shearing force acting 


FIG. 3.—Analysis Model for Multilayer Diaphragm System 


between layers, and and and the extreme fiber distances from 
the centroids layers, and respectively. 

The preceding approach two-layer systems can extended include 
multilayer systems. Referring Fig. the slip between connected layers gives 
rise the forces shown. The strain the bottom layer, given 


which the elastic modulus; A,, cross-sectional area; and moment 
Similarly the strain the top fibers layer, can written as: 


Fy Fix 
FIG. Model for Two Layer Diaphragm System 
3 Fj i Cj Fij i 
Fim F c F Fim 
Fim m m Fim Fmn 
Fin n Ca Fn 
by: 


The rate change slip between layers, and along the length the 
diaphragm given the difference strains Eqs. and Thus: 


which the stiffness glue line which function the modulus 
rigidity, thickness, and width glue. 
Upon substituting Eqs. and into Eq. the following obtained: 


M,C, M,C, 


Assuming that each layer bent the same curvature and upon substituting 
the values and terms all couples into Eq. expression 
obtained relating the rate change slip along the length the diaphragm 


and the shear forces. This procedure repeated for each glue line obtain 
the following: 


and each layer; and the number layers the diaphragm. 

Eq. contains systems second-order linear differential equations that must 
solved simultaneously any given point the beam determine shear 
forces between adjacent layers. Once shear forces are determined, the moments 
and stresses each the layers are obtained. 


AND 


Because the expense associated with the construction and instrumentation 
lumber diaphragm, only one was tested. The test was carefully conducted 
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and measurements various component properties were obtained. 

The test diaphragm (Fig. was (4.57 long and (1.52 wide. 
consisted layers that were glued with Scotch-Grip Wood Adhesive 
No. 5230. The size each layer was 1/2 in. 1.5 in. (12.7 38.3 
mm). 

Eleven 7/8-in. 3-1/4-in. (22.2-mm 82.5-mm) joists were used spacings 
1.5 (0.46 center. These joists provided the overall lateral stability 
the diaphragm. Each joist was approx (1.52 long. 


- 


The 15-ft long lumber layers were installed starting from one end the 
joists. Each layer was nailed its center each joist. Before placing the 
next layer, nails approx 0.09-in. (2.3-mm) diam were placed act glue 
line spacing shims keep the glue line thickness constant (Fig. 5). Glue was 
then applied using cartridge applicator and the next layer was placed and 
nailed the joists. All the layers were placed the same manner. After the 
glue was cured, any excess glue outside the plane the diaphragm was taken 
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off using chisel. The final dimensions the diaphragm are: (1) Length 
(4.5 m); (2) width (1.5 m); (3) average thickness glue line 
0.09 in. (2.3 mm); (4) average width glue line 0.5 in. (12.7 mm); (5) 
average width lumber layer 1.9 in. (48.25 mm); and (6) average thickness 
lumber layer 0.48 in. (12.2 mm). 

The wood used for decking and joists was western hemlock. Moisture content 
was about 12% average. static load test method was used measure the 
modulus elasticity each wood strip before diaphragm assembly. Each 
piece wood was supported each end with known span length. midspan 
small known load was placed and the deflection due that load was measured 
with dial gage. Knowing the span, support conditions, load, deflection and 
moment inertia the section, the modulus elasticity, was calculated 
for each piece wood. 

elastomeric adhesive, Scotch-Grip Wood Adhesive No. 5230, was used 
for bonding the diaphragm resist shear between the decking boards and thus 
increase the diaphragm stiffness. The adhesive contains solvent and cures 
without much shrinking when exposed atmospheric moisture. The resulting 
seal water resistant and flexible, yet retains considerable strength and stiffness. 

The shear modulus the adhesive important for this study and was measured 
testing samples double shear. 

The diaphragm was placed frame (Fig. 4). was simply supported 
both ends. Some guide beams were provided the test frame prevent lateral 
deflection buckling shown Fig. 

uniformly distributed load 200 was applied system 
weights and pulleys. Eight pulleys were used which were equally spaced 
and transferred the gravity load the diaphragm (Fig. 4). These small equal 
concentrated loads may considered uniformly distributed load with 
negligible error. 

Before taking any reading, the diaphragm was loaded 200 and unloaded 
times. The reason that shear modulus glue decreases cycles 
loading (6). Loads were applied the diaphragm steps. Increments 
load were chosen arbitrarily. After each increment was reached, 
readings were taken and the load was released zero again. The load was 
then increased the next greater load. 

Deflection was measured scale with 1/32-in graduations. The scale was 
fixed the middle section the diaphragm with mirror its side. taut 
cord was stretched from the two ends the diaphragm passing across the 
scale and the mirror. When the diaphragm was not loaded, the reading the 
cord the scale was taken. the diaphragm was loaded, the scale moves 
with the diaphragm but the cord remains its position. The scale read 
each load and the difference two readings gives the deflection the diaphragm 
for that load. The mirror helps prevent any parallax reading. 

Strains different layers the diaphragm were measured strain gages. 
Strain gages were placed the middle the layer measure the average 
strain. Because the strain the top and bottom fibers any layer are not 
the same, strain the middle that layer gives the average axial strain. Strains 
alternate layers were measured. Positions strain gages across the diaphragms 
are shown Fig. 


analyze the diaphragm analytically, the slip modulus, defined 
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Gb/t, needed. Here, the shear modulus glue, and are the 
width and thickness glue line, respectively. The actual slip modulus glue 
lines the diaphragm were measured taking double shear specimens from 
the diaphragm. The diaphragm was cut and the specimens were taken after 
all other tests were completed. 


The following sequence tests was conducted: 


The first test was conducted for comparison with the results obtained 
from the theoretical analysis Ref. the first test, the diaphragm had 


Diaphragm 


Strain Gauge No. 


FIG. 6.—Positions Strain Gages Diaphragm 


boundary members discontinuities the layers. Each decking layer was 
connected each joist single nail that resisting moment could 
develop. This accordance with the theoretical method presented Ref. 
The diaphragm was loaded 200 and the readings deflection 
and axial strain were recorded. 

The second test measured the realistic effect discontinuous layers. The 
same diaphragm was used for the second test. The layers the diaphragm 
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were cut transversely with sabre saw, produce end joint pattern typical 
real diaphragms (Fig. 7). The lengths different pieces the courses 
lumber decking were selected using the limitations selecting random lengths 
described the American Institute Timber Construction (AITC) Standard 
112-65 found Ref. 14. The diaphragm was loaded and readings deflection 
and strains were made. 

The third test was done study the effect gluing the joists the 
deck courses the deflection and strains the diaphragm. The same diaphragm 
was used for the third test. The joists were glued with the same glue. Thirty 
days were allowed for curing the glue. Following the same loading sequence 
the deflection the diaphragm and the strains different layers were measured. 


FIG. 7.—Discontinuities Decking Diaphragm 


The fourth test was designed study the effect boundary members 
the diaphragm. the four boundaries the diaphragm four boundary 
members were glued. Thirty days were allowed for glue cure before testing. 


addition the analytical results, four different sets experimental results 
were obtained from this study. All these results consist strains different 
layers the diaphragm and deflection midspan. 

The first set consists results obtained from the theoretical analysis 
the diaphragm, using the computer program DAD—Difference Analysis 
Diaphragms (5). The computer program does not account for any discontinuity 
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the decking the diaphragm and does not account for any effect from 
the joists boundary members. 

The second set results was obtained from the test the diaphragm model. 
The diaphragm for this test had discontinuity its decking and boundary 
members. The joists the diaphragm were not structurally connected the 


TABLE 1.—Deflection Obtained from Theoretical Analysis and Model Test 


Variable 
(1) 


Theoretical analysis 
Model test 


Note: Loading 200 uniformly distributed; 4.448 in. 25.4 mm. 


Deflection mid-span inches 


(2) 


TABLE 2.—Strains Different Layers Obtained from Theoretical Analysis and Model 
Test 


Layer number 


diaphragm Theoreticai 
starting from top analysis Model test 
Gage (total number strain, micro strain, micro 
number layers 31) inch per inch inch per inch 


(1) (2) (3) (4) 


“No reading obtained. 


Note: Positive strain indicates tension; Loading 200 length; 4.448 
in. 25.4 mm. 


diaphragm (one nail was driven the intersection each joist and deck course). 

This model had identical properties the diaphragm analyzed theoretically. 
The third set results was obtained testing the model with discontinuities 

its layers. The comparison between the second and third set results shows 

the effect discontinuities the layers diaphragm. 

The fourth set results was obtained testing the model with discontinuities 


0.735 
—330.7 —288.23 
—9.50 
+2.50 
+0.1 +0.0 
+0.1 +7.33 
+0.1 +21.17 
+14.5 +27.17 
+67.9 +64.50 
+334.8 293.33 
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its layers and with the joists glued the decking, thus developing resisting 
couples. 

The fifth set results was obtained testing the diaphragm with boundary 
members all its four sides. Boundary members were 0.5 in. 1.9 in. 
cross section. 

The comparison the fourth and fifth sets with the third set shows the 
effects joists and houndary members the diaphragm behavior. 

Tables and show the first and second result sets. These tables indicate 
that theoretical and experiment results are reasonable agreement. The theoreti- 


200 


Experimental 
Theoretical 


COMPRESSION TENSION 


FIG. 8.—Strain Distribution Diaphragm 


cal analysis overestimates deflections and axial strains the critical outer layers. 
The deflection obtained the theoretical analysis about 10% more than 
the experimental value. Outer layer strains obtained theoretical analysis are 
14% greater than experimental strains. Strain distribution pattern was almost 
identical both cases (Fig. 8). 

The reason for this discrepancy that the theory, its development, does 
not account for the effect joists. The joists were nailed the decking. 
was assumed that the joists act like stiffeners plate girder and distribute 
load the layers well prevent buckling. However, for lumber diaphragm 
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Type diaphragm 

(1) 

Diaphragm without any discontinuity its 
layers 

Diaphragm having discontinuous layers 

Diaphragm having joists glued 


members 


Conditions 


with 


(2) (3) 


—288.23 


+2.50 
+0.0 
+7.33 
+21.17 
+8.50 
+5.67 
+27.17 
+64.5 
+293.33 


(bottom) 


“No reading obtained. 


in. 25.4 mm. 


Diaphragm having glued joists and boundary 


+350.17 
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where the courses decking slip with respect each other, the joists develop 
some resisting moment. other words, the joists oppose the interlayer slip 
the decking and make stiffer than the theory assumes. result the 


Diaphgram 


with 


continuous discontinuous 


layers layers 
Layer strain, strain, 
number micro inch 
from top per inch per inch 


(4) 
—322.14 
—21.33 
—2.50 
+5.0 
+6.83 
+18.33 


+13.33 
+31.83 
+65.83 


TABLE 3.—Deflections Obtained from Model Tests for Various Conditions 


Deflection midspan, inches 


(2) 


TABLE 4.—Strains Different Layers, Obtained from Model Tests for Various 


Diaphgram 


with 


glued joists 
strain, 

micro inch 

per inch 


(5) 


+11.33 

+26.5 

+59.0 
+339.67 


joists reduce the deflection the diaphragm. They also affect the strain 
distribution, reducing the outer layer strains and increasing the inner layer strains. 

The average modulus elasticity layer was measured. actuality the 
modulus elasticity wood layer varies along its length. The experimental 
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Diaphgram 
with bound- 
ary mem- 
bers glued 
strain, 
micro inch 
per inch 


(6) 


Note: Positive strain indicated tension; Loading 200 lb/ft length; 4.448 


0.73 
0.78 
0.64 
0.63 
number 
(1) 
—0.83 
—2.17 
—16.83 
+14.25 
+20.0 
+163.0 
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strains are local strains while theoretical strains are average values. This 
another source error which can further evaluated with additional tests. 

Tables and show the second, third, fourth and fifth set results. These 
tables compare the effect discontinuities the layers and the effect joists 
and boundary members. They also show the deflection the diaphragm increases 
about and the outer layer stresses increase about 15% due the discontinuity 
its layers. 

The effect discontinuities not very significant. When the joists were 
glued the decking, the deflection the diaphragm was reduced about 18% 
(Table 3). The change strain was insignificant (Table 4). These are important 
observations for the designer, and were less than expected. 

The effects boundary members the deflection the diaphragm was 
insignificant but reduced the outer layer strain about 50% (Tables 4). 

The reduction the outer layer strains due the fact that boundary 
members become the outer layers the diaphragm and change the outer layers 
inner layers. 

There significant stiffness gain from the addition boundary layers 
and should possible compute their effect deflections incorporating 
them into the analysis additional layers. The small effect the boundary 
members may seem unusual since they are usually considered very important 
conventional plywood sheathed diaphragm design. Their small effect due 
the fact that the boundary members were not large, especially when compared 
the size the lumber deck members. This special feature thick 
lumber diaphragms. plywood sheathed diaphragms, the relatively large bound- 
ary members contribute more the overall performance. 


Conctusions 


The experiment verified the theoretical analysis diaphragm satisfactory 
manner. The study also defined the effects different structural features 
the diaphragm. From this study the following conclusions may drawn: 


The theoretical analysis the diaphragm over estimates strains and 
deflections because does not account for the effects joists the diaphragm. 

Presence discontinuities layers glued lumber diaphragm does 
not have considerable significance. 

Gluing the joists the decking lumber diaphragm reduces deflections, 
but has little effect strains each layer. 

Addition boundary members does not have significant effect the 
deflection parallel course lumber diaphragm for members the size used 
this study. Larger boundary members might have greater effect, but probably 
boundary members contribute less diaphragm performance for glued lumber 
diaphragms than would the case for nailed plywood diaphragms. 

Additional testing diaphragms needed reach more conclusive results. 
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The following symbols are used this paper: 


cross section area layer 

distance from centroidal axis layer extreme fibers that layer 

force between layers and 

coefficient matrix; 

external moment; 

total number layers; 

shear connected stiffness between layer and 
deflection; 

strain; and 

stress. 
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EFFECT BOWING RECTANGULAR 
PLANE FRAMES 


When member subjected bending, there shortening the chord 
length, shown Fig. This change length has been called bowing 
deformation curvature shortening. Researchers the area inelastic frame 
analysis continue debate the merits including bowing second-order analyses 
(1). There are structural configurations for which thorough nonlinear analysis, 
including bowing effects, necessary describe structural response. But, 
bowing small comparison axial deformation that the effect negligible 
changing rectangular frame response? 

The purpose this paper present results frame analyses that allow 
the reader compare the magnitude bowing and axial deformation, and 
examine the influence bowing altering structural response service 
and collapse load level (4). The analyses are based second-order elastic-plastic 
analysis (beam-column rotational and translational stiffness with allowance for 
the formation plastic hinges). addition, axial deformation and bowing are 
included (4,5). 


The change chord length, 5,, caused bending shown Fig. Bowing 
deformation evaluated the following formula (3,7): 


The mathematical expressions for have been presented several authors 
(e.g., see Ref. 6). The expressions used herein are shown Appendix 
function end moments, and axial load, original member length, 
and bending stiffness, 


For the case central lateral load, bowing evaluated parts (see 
Fig. 2). The reduced chord lengths, and are computed separate 
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bowing analyses. Denoting the angle inclination the subchords the 
main chord the overall bowing deformation given 


The aforementioned formulation also valid for the case plastic hinge 
either end the load point. 


FIG. 2.—Bowing Parts 


Note that bowing deformation dependent the squares member slopes. 
Accordingly, the use small deflection assumptions for seem inappropriate 
analysis considering bowing. However, when bowing deformation, 
compared deformation caused axial load, the relative magnitude 
terms can only judged examining typical frame analyses. 


Frame 


Several single-bay, multistory frames have been analyzed collapse (both 
with and without the consideration bowing) (4). Thus, possible compare 
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the magnitude bowing deformations with axial deformations service and 
collapse load level. Also, possible examine how much the collapse capacity 
and overall frame deformations are influenced bowing. 

The results all the analyses are well represented the behavior 
the four-story frame shown Fig. The frame strong-column”’ 
design suggested Horne and Majid (2). The frame was proportionally loaded 
collapse. Thus, any load level represented load factor, ranging 
from unity service loads the collapse load factor. The complete history 
load, deformation, and plastic hinging was recorded (4). The summary 
member shortenings (5,, axial, and bowing) are shown nondimensional 
form Table 

Behavior Service Loads, 1.—For the imposed loading, the columns 
tend bent into double curvature, while transversely loaded girders have 
much higher contribution single curvature. With reference Eq. Appendix 


F/2 


GIRDERS: W16X40 
COLUMNS: 
W12X79 (BOTTOM STORY) 
W10X60 (ELSEWHERE ) 
12' (3. 66m) 
30'(9.15m) 
3.22t(32.1 KN) 
13.4t(133.6 KN) 


TON 2.24 KIPS) 


FIG. 3.—Frame, Properties, Service Loads 


can shown that for equal forces the ratio bowing deformation for 
double curvature never more than fifteenth the bowing for single curvature. 
Thus, for the girders should higher than for the columns. Also, 
known that the columns have the larger axial strains. Both observations 
were verified the results: (1) The axial deformations the columns was 
generally two orders magnitude higher than the bowing the columns; and 
(2) the bowing and axial deformations the girders were the same order 
magnitude. Even when both bowing and axial deformations are considered, 
the total nondimensional shortening any girder was less than that for the 
least compressed column. Realizing that axial deformation secondary 
bending, the bowing effects service load are clearly negligible. 

Behavior Collapse Loads, 1.499.—The order plastic hinge formation 
shown Fig. Failure occurred instability prior the formation 
complete collapse mechanism. Each girder had developed plastic hinges 


|| F F/2 
H/2 
Fr2 F 
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TABLE 1.—Member Shortenings 


Left- Column Right-Side Column 


Top 


Bottom 


Top 
Bottom 
Note: axial deformation bowing deformation (see Appendix 
are positive when chord length shortened. 


FIG. 4.—Frame Collapse Load Level (Numbers Indicate Order Hinge Formation) 


(1) (4) (5) (6) (7) 
(a) Service Loads, 
0.0559 0.0004 0.0420 0.0100 0.0608 0.0012 
0.1081 0.0002 —0.0036 0.0065 0.1252 0.0006 
0.1557 0.0000 0.0119 0.0083 0.1942 0.0010 
0.1541 0.0013 0.0066 0.0080 0.2006 0.0007 
(b) Collapse Loads, 1.499 
4 i 
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the center and right end, and earlier formed hinges the bottom two stories 
had time develop appreciable rotations. Now, bowing terms predominate 
for the girder, and are the same order magnitude the larger axial 
shortenings the columns. Bowing deformations the columns were negligible. 

The only time that bowing deformation was appreciable was when plastic 
hinges developed midspan, and had undergone rotation. Yet, the result 
such deformation the overall frame behavior was negligible collapse! The 
collapse load factor, changed only 0.001 (from when bowing 
was considered. The top sway height ratio, A/h, was 0.0160. When bowing 
deformations were also included, the ratio changed 0.0001. 


The results obtained for the four-story frame were typical for all frames 
analyzed. The state deformation collapse still sufficiently small, 
that the basic mode frame behavior predicted simplified models adequate. 
This computer experiment, though limited scope, does support the contention 
that bowing deformations are negligible rectangular planar frames. 


The bowing deformation, 5,, shown Fig. has been evaluated the 
following formulas (4): 


where and member length and bending stiffness, respectively. The 
bending moments, and are positive they are counterclockwise. The 
stability factor, function axial load, and always positive. The 
aforementioned expressions are for compressive axial load. the case tension, 
hyperbolic functions replace the trigonometric functions. 

For small axial load 0.5), the following series expressions are 


inwhich 
26°/45 


360 


Where dual signs are shown, the positive sign applies for tensile axial load 
and the negative sign for compressive axial load. 
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DISCUSSION 


The writer thanks Jain and Hays for their enlightening discussion. Jain’s 
observations are answered here first. The subject matter the paper forms 
only part the writer’s more general work (17), and the strut inelastic behavior 
one the examples chosen for presentation. Thus not all the literature 
strut modeling has been included. few the papers mentioned Jain are 
quoted Ref. 17. few others are almost contemporary period. 

The convergence difficulty mentioned the paper was specific the type 
discrete element (6,17), and not due formation plastic hinge midspan. 
the interior each discrete element plastification monitored two cross 
sections. Thus when the kink positioned exactly midspan, two plastic 
hinges are formed, due symmetry, close the midspan, one each 
the two discrete elements either side. 

This can overcome defining behavior for one 
these two discrete elements. Such procedure has been used Harung 
and Millar (42) overcome mathematical singularity resulting from the 
simultaneous formation two plastic hinges structural node where two 
members with equal plastic properties meet. This will, some extent, affect 
the accuracy the predicted deformations and therefore was not adopted. 
Adoption the noncentral kink was used better contrivance. There are 
the other class discrete elements where plastification occurs the junction 
two elements. such models, the kink can positioned midspan. The 
writer wishes add that considerable number days were spent trying 
debug, before the aforementioned phenomenon two simultaneous hinges 
was spotted the cause convergence failure. Refs. and contain the 
remedial measures that can adopted overcome convergence difficulties 
due inelastic material behavior. Ref. contains example convergence 
failure strut inelastic analysis due the combined effect geometric 
nonlinearity and the method discretization. The problem here the kink 
originally input simulate imperfection removed due inelastic reversals. 
The subsequent analyses are thus made closer the difficult problem buckling 
perfect structures. 

Jain quotes Kaldjian and Fan (30) the closeness the results obtained 
Ramberg-Osgood and elasto-plastic moment-curvature models. This not 
surprising. The Ramberg-Osgood moment-curvature relation for cross section 
was derived from Ramberg-Osgood representation the stress-strain relation. 
This accounts only for the gradual plastification across the depth the section, 
against the abrupt transition inherent the ideally elastic-perfectly plastic 

1979, Santhanam (Proc. Paper 14333). 
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moment-curvature relation. The more important effect the strain-hardening 
mild steel not accounted such procedure unless the Ramberg-Osgood 
algorithm modified and enhanced, for example, done Eisenberg (41). 
The normal Ramberg-Osgood model and the multilinear representation the 
Masing model (7,38) may said equivalent. The advantage the latter 
that other types stress-strain, general, force-deformation behaviors, 
such shown Table for example, can easily incorporated. The improved 
model Eisenberg encompasses wider area including multi-axial states 
stress. 

The writer does not wholly agree with Jain’s conclusion the proposed 
stress-strain model. Inelastic modeling techniques are necessarily computationally 
expensive, and thus, are not directly adopted for analysing large structural 
systems. Some the alternatives adopted dealing with such structures are 
restrict the material model sophistications small region and use finer 
mesh size, adopt condensation techniques separate regions requiring more 
rigorous treatment. either these good apriori knowledge the location 
the probable influential zone the structure needed. dynamic analyses, 
not static analyses, such guesses are difficult make, and wrong judgement 
can lead significant wasted effort. The writer proposes third approach 
wherein lower-order macro-analysis the entire structure performed first, 
and subsequent more sophisticated micro analysis one two selected 
smaller segments are carried out for further detailed study, e.g., ductility 
requirements, etc. The boundary conditions, including their history, the cutoff 
nodes defining the smaller segment can obtained from the historical displace- 
ments, accelerations etc. the macro-analysis stored tape. This treatment 
differs from the usually understood substructure analysis that based 
hindsight and not dependent foresight. 

program such FRAME can successfully used analyze moderately 
sized structural segments about 20-25 joints and 15-20 inelastic members, 
though the total number members can more (43). Noting the organization 
the computations wherein the member inelastic analyses are isolated from 
the global analysis the structural nodes (6,17), structural segment about 
total 3,000 degrees freedom-4,000 degrees freedom can realistically 
solved. The band width greatly reduced compared that frontal solution 
similar structure. Schrem (44) considers such structure decent 
size for inelastic analysis though not for linear elastic analysis. The writer 
has found that the simpler elasto-plastic model predicts very large ductility 
demands the local regions for which might difficult design the structural 
connections. The writer agrees with the general spirit behind Jain’s conclusion 
that attempts must made simplify material models the light more 
experimental data, consistent with cost-accuracy considerations. 

The writer thanks Hays for comparing the present model with Masing type 
treatment (7,38) the virgin stress-strain curve mild steel. This Masing 
representation was one the early models tried the writer towards his 
research (17) but was discarded because had correlation with experimentally 
observed behavior the reversed range. Such model was used Chen 
and Atsuta (40) for stress-strain representation wherein gross under estimation 
could seen the negative side the stress axis. The same behavior has 
been reported Hays whose model bridged about 18% gap between the 
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experiment and the model the final cycle Fig. but undershot 
much 36% the negative side stress (39). Thus the adoption Masing 
decomposition for the special shape the virgin stress-strain curve mild 
steel yields good fit for monotomic loading, despite one the three basic 
components being negative stiffness (39), which difficult give reasonable 
physical meaning to. However this deficiency persists under reversed conditions 
quite contrary experimentally observed behavior (14). This explains the 36% 
deviation from the test results and the inaccuracies the locally predicted 
moment-curvature responses (39). 

Accepting these inaccuracies and deviations, computationally more cheap 
adopt two-component Masing curve the shape shown Fig. and 
analyze the structure Fig. than use Hays’ three-component Masing 
curve. This avoids the convergence failure represented Fig. 11, but overshoots 
the experimental curve the positive side the load axis, can easily 
expected. The writer grateful Hays for his observation that the proposed 
a-B model would give more precise results for other load histories. Hays has 
also identified probable increase the base rotational stiffness the frame 
Fig. one the reasons for the increased stiffness and load Fig. 
the reversed direction (39). 


40. Chen, F., and Atsuta, T., Response Column Segments under Biaxial 
Journal the Engineering Mechanics Division, ASCE, Vol. 99, No. EM4, 
Proc. Paper 9957, Aug., 1973, pp. 685-701. 

41. Eisenberg, A., Generalization Plastic Flow Theory with Application 
Cyclic Hardening and Softening Journal Engineering Materials and 
Technology, Vol. 98, July, 1976, pp. 221-228. 

Harung, S., and Millar, A., Failure Analysis Skeletal Plane 
Journal the Structural Division, ASCE, Vol. 99, No. ST6, Proc. Paper 9776, June, 
1973, pp. 1051-1074. 

Santhanam, K., Dissertation Abstracts International, Vol. No. 1978. 

Schrem, E., Aspects Software Systems for Nonlinear Finite Element 
Finite Elements Nonlinear Mechanics, TAPIR, Norwegian Institute 
Technology, Trondheim, Norway, 1978. 


Errata.—The following corrections should made the original paper: 


Page 199, paragraph line Should read mild steel 

Page 202, paragraph line Should read ‘‘The spotted line represents the 
instead dashed line represents the 

Page 202, Fig. Line segments and must shown equal 
slope and not printed. 

Page 217, Appendix central block the flow chart: Should read o,,, 
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PARAMETRIC INVESTIGATION CABLE 


The authors presented theoretical investigation vibrating cable networks. 

scale factor needed comparing model with the real structure. 
the material the model the same that the real structure, the cross-sectional 
areas the cables the model must such that the stresses the cables 
both the model and the prototype are identical. Thus, the statement 
while leaving the cable areas mentioned the authors 
471, unrealistic. For real suspension structure, not logical build 
laboratory model using the same cable areas the real structure. Also, the 
use the same cable areas (and, course, the same sag-span ratios) will 
mean that the load intensity the model must greater than that its 
prototype order obtain same cable stresses (assuming same cable material). 
But, again, this not logical. the load intensity (per unit area the surface) 
kept unchanged, the cable areas must altered according the relation 


which cross-sectional area the new system (model); and 


cross-sectional area the base system 
Also, Eq. incorrect and should have the form 


This true, again, provided that the intensity the total (dead live) 
loads per unit area the surface identical both the model and the real 
structure. 

Eq. may correct for only single cable (not network), the load 
intensity per unit length the same both the model and its prototype. However, 
for network, Eq. may correct 


which load intensity per unit area the model; and load intensity 
per unit area the prototype. This not logical. more reasonable choice 
may 


This gives network 
The use 18, 20, and gives 


“March, 1979, Louis Geschwindner, Jr. and Harry West, (Proc. Paper 14436). 
Faculty Engrg., Ain Shams Univ., Abbassia, Cairo, Egypt. 
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while the use Eqs. 18, 20, and yields 


Both 43a and are obviously different from Eq. 23. However, the 
authors did not offer any applications that can prove the correctness 
and 23, and Figs. and represent the effect the corresponding parameters 
the dynamic response cable network rather than single cables. 

472, the authors said: 


complete the investigation the influence applied load the 
natural frequencies and mode shapes the network, equal variation 
force applied each joint the network without the corresponding 
change being made the mass the system. 


this statement, the writer has the following comments: 


The natural frequencies and mode shapes are obtained for free vibrations. 
Thus, exciting forces are needed. 

The application any dead live loads joints result only change 
the masses these joints, but not the forces (which are zero for free 
vibrations). 

The fourth conclusion given the authors 476, well the 


work Krishna, al. (8) support the point view the writer mentioned 
his second comment. 


should noticed that, due the nonlinearity cable structures, the 
stiffness (or flexibility) matrix different for any individual loading condition 
than for any other loading condition. secant stiffness (or flexibility) matrix 
for the working loads condition gives usually reasonable results (20). 

his work, the writer (20) uses the flexibility matrix instead the stiffness 


which the flexibility matrix [K]~' obtained directly and not finding 
the inverse the stiffness matrix [K]. 

Finally, for Fig. 4(a) and the corresponding section pp. 472-473, the word 
must replace the word The change loading yields corre- 
sponding change both the stiffness and mass matrices the cable network. 
seems the writer that the authors were confused between the exciting 
forces and live loads. Exciting forces produce forced vibrations, while live 
loads are added the dead loads the mass matrix. 


20. Hegab, and Dynamics Suspension Bridges with Notes Suspension 
Cables and thesis presented the University Tasmania, Hobart, Australia, 
1977, partial fulfillment the requirements for the degree Doctor Philosophy. 
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Continuous 


The writers would like thank Starkey and Rabb for their comments regarding 
the role they feel that limit analysis should have the design bridge structures. 
The writers agree that designing each section with equal overload capacity, 
the importance longitudinal redistribution moment continuous span 
structures will diminished. However, the ultimate capacity each section 
depends the post-elastic effective width the slab section. Also, main 
conclusion the paper—namely, that even prior yielding, the distribution 
factors the positive and negative moment regions differ—will have effect 
the ability design create simultaneous yielding each section 
for the same overload. Despite these questions, the writers concur that, 
present, the design approach outlined Starkey and Rabb proper. 

With respect the comments and questions raised Kostem and Hall, 
the writers have reviewed the experimental data for the three-span model 
order clarify the effect support movement the reported displacements. 

The support movements recorded during testing are given Table 

the paper the reported displacements the specimen the center line 
the loaded span, were computed follows: 


which measured midspan displacement; and average measured 
displacement interior supports. The use Eq. does not account for the 
uplift the exterior supports. moment distribution analysis indicates that 
the following relation more accurately reflects this influence: 


which average measured uplift exterior supports. 

Using Eqs. and and the data given Table yields the corrected 
displacements presented Table which indicates that consideration the 
exterior support uplift reduces the estimated bending displacements. 

The main purpose the paper was interpret the test results they relate 
the distribution load and the effective widths the slab. These factors 
were investigated using strain data and revealed that the effective widths were 
in. (0.2 and in. (0.15 for the center girder loads kips (90 
and kips (180 kN), respectively. Corresponding values for the edge 
girders were 19.4 in. (0.5 and 16.5 in. (0.4 m), respectively. The load carried 
the center girder was 0.55 kips), and 0.58 kips). 

“September, 1979, Paul Botzler and James Colville (Proc. Paper 14824). 

Struct. Engr., Tamarack Corp., Glen Arm, Md.; formerly Grad. Research Asst., Dept. 


Civ. Engrg., Univ. Maryland, College Park, Md. 20742. 
Dept. Civ. Engrg., Univ. Maryland, College Park, Md. 20742. 
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Using these results, simple elastic theory gives the values given Col. 
Table These latter values compare well with the adjusted displacements 
computed using Eq. 


Results reported the discussers for kips (90 kN) are: 


BOVA—center girder: 0.0343 in. (0.87 mm); edge girders: 0.0132 
in. (0.34 mm). 


SAP4—center girder: 0.0304 in. (0.77 mm); edge girders: 0.0126 
in. (0.32 mm). 


These results suggest that the load carried the center girder around 
0.547 P-0.565 This compares favorably with the value 0.55 given previously 
the writers. Furthermore, using effective width slab in. (0.53 
m), and load 0.55 the center girder, simple theory indicates that 
load kips (90 kN), the center girder would displace 0.0314 in. (0.80 


TABLE 5.—Three-Span Model: Support Movements, inches 


Load, kips 
Support location 


(1) 


(a) Exterior supports (average) 


Center girder —0.007 
Edge girders —0.004 


Interior supports 


Center girder (left) 
Center girder (right) 
Edge girder (left) 
Edge girder (right) 
Edge girder (left) 
Edge girder (right) 


Note: kip 4.45 kN; in. 25.4 mm. negative displacement indicates upward 
movement; positive displacement indicates downward movement. 


TABLE 6.—Three-Span Model Bending Displacements, inches 


Simple elastic 
Location Load, kips theory 
(1) (5) 
Center girder 
Edge girder 
Center girder 
Edge girder 


“Displacements computed using simple elastic beam theory and effective widths and 
distribution factors from experimental strain data. 
Note: kip 4.45 kN; in. 25.4 mm. 


(3) 
+0.013 +0.023 
+0.008 +0.015 
+0.012 +0.019 
+0.014 +0.022 
+0.009 +0.012 
+0.011 +0.025 
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mm) and the edge girder displacements would 0.013 in. (0.33 mm). These 
values are close agreement with the numerical results reported Kostem 
and Hall. is, therefore, suggested that the major difference between the 
numerical results obtained the discussers and the experimental data presented 
the writers not due the effect the support movements, but rather 
due the reduced effective width slab acting compositely with the steel 
beam. 

literature survey Fan (28) related effective width indicates the 
complexity the problem, producing corresponding wide variety equations 
and procedures for predicting effective widths for use with simple beam theory. 
brief review these procedures indicates that predictions the effective 
width slab the three-span tests reported the writers could range from 
in. (0.15 m)-18 in. (0.46 m), and that stress concentrations due concentrated 
loads will result smaller effective widths than those induced uniform 
loads. These indications are compatible with the findings the writers. 

summary, since has not been possible for the writers perform 
review the numerous assumptions made Kostem and Hall their numerical 
analysis procedures, the writers are unable explain satisfactorily the reasons 
for the reported discrepancies between test and theory. However, the writers 
believe that unlikely that these differences are due solely the effect 
support movements during testing. more thorough consideration these 
effects reduces but does not eliminate the noted differences. suggested 
that the relatively small width concrete acting compositely with the steel 
beam may fact the primary reason why the experimental displacements 
are noticeably larger than those predicted the numerical procedures. 

The writers would pleased discuss this possibility further with Kostem 
and Hall. 


28. Fan, Width Composite Bridges Ultimate thesis presented 
the University Maryland, College Park, Md., 1974, partial fulfillment 
the requirements for the degree Doctor Philosophy. 


DISCUSSION 


Closure Egor ASCE, Victor Zayas,° 
and Stephen Mahin,’ ASCE 


The writers would like thank Hays for his interest the paper and for 
raising important question regarding material properties the model tubes 
they relate actual structures. This question has been carefully considered, 
and can seen from Fig. the paper, experimental results with unannealed 
tubes are totally misleading. the other hand, the stress-strain diagrams for 
annealed specimens closely resemble those for steels used real structures. 

The writers believe that fully annealed specimens provide reasonably accurate 
models for simulating the inelastic cyclic behavior tubular members used 
the field. would impractical and prohibitively expensive fabricate 
4-in. tubing using the same techniques are used manufacture 24-in. 
rolled pipes obtain reasonable simulation residual stresses. 
recognized that initial residual stresses will affect the first buckling load 
column. These effects, however, can accounted for the initial buckling 
load predicted using the stress-strain curve from column stub test 
full size member (16). The effects the initial residual stresses the buckling 
load the second and subsequent cycles usually greatly reduced because 
compressive and tensile axial yielding during the first inelastic cycle. The 
writers recommend that for the inelastic analysis braced frames, data for 
both monotonic and cyclic material properties the members the frame 
used predict the initial and cyclic buckling loads (17). 


16. Guide Stability Design Criteria for Metal Structures, Johnston, ed., 3rd 
ed., John Wiley and Sons, Inc., New York, N.Y., 1975. 
17. Zayas, A., Popov, P., and Mahin, A., Inelastic Buckling Tubular 


Steel Report University California, Berkeley, Calif., 
June, 1980. 


1979, Egor Popov, Victor Zayas, and Stephen Mahin (Proc. 
Paper 14982). 

Civ. Engrg., Coll. Engrg., Div. Struct. Engrg. and Struct. Mechanics, 
Univ. California, Berkeley, Calif. 94720. 

Student Civ. Engrg., Univ. California, Berkeley, Calif. 94720. 

Asst. Prof. Civ. Engrg., Univ. California, Berkeley, Calif. 94720. 


MARCH 1981 


PracticaL Cope 


The writers would like thank Carskaddan and Haaijer for their discussion 
one the aspects the OHBD Code calibration. 

the OHBD Code, kinking defined permanent deflection, and the 
tolerance set 0.001 the span length. The code allows load causing 
permanent deformation occur average once yrs. Load factors were 
derived accordingly. 

Several steel girder bridges, recently built Ontario, were considered. The 
allowable permanent deflection (kinking) never governed the design. seems 
that the differences between Carskadden and Haaijer’s results and the writer’s 
findings are due different definitions kinking and tolerancies. 


The writers would like thank Tabba for his discussion and contribution. 
response the particular questions Tabba has noted, offer the following: 


thorough study (12) typical bridges, constructed the United States, 

indicates that the central angle between supports always less than 40°. Thus 

the condition 45° not practical situation and thus not relevant. 
The definition low torsional stiffness has not been defined Tabba; 

however the bridges under study have significant torsional stiffness, 

indicated Table 
The bridges under study have symmetrical cross sections, thus coupling 

between flexural vertical and torsional vibrations valid. 


1979, Andrzej Nowak and Niels Lind (Proc. Paper 15061). 

Prof., Coll. Engrg., Dept. Civ. Engrg., Univ. Michigan, Ann Arbor, 
Mich. 48109. 

Dept. Civ. Engrg., Univ. Waterloo, Waterloo, Ontario, Canada. 

December, 1979, Conrad Heins and Sahin (Proc. Paper 15064). 


“Prof., Dept. Civ. Engrg., Inst. for Physical Science and Tech., Univ. Maryland, 
College Park, Md. 20742. 


Student, Dept. Civ. Engrg., Univ. Maryland, College Park, Md. 20742. 


ST3 DISCUSSION 587 


The example relating the computed frequencies for two-span bridge relative 
single-span bridge were not computed correctly. 


The results should 1.006 0.62 and 0.43 which conforms 
the general trend. The discusser should have used the referenced maximum 
span length (NL) 21. all cases the maximum span length should 
used applying Eq. 21. 


12. Heins, P., Girder Bridge Design—State Engineering Journal, 
American Institute Steel Construction, Vol. 15, No. Dec., 1978. 


Current STRUCTURAL ENGINEERING REFERENCE 
Closure Sam 


The writer would like thank Saul for having pointed out the discrepancies 
the reference sources listed the paper and submitted corrections them. 
especially appreciates the efforts the reviewers his paper because 
they have done tremendous job the process. Had not been for their kind 
efficiency, the writer would have still been revising the manuscript instead 
writing closure like this. Last, but not least, the writer’s sincere 
gratitude toward the ASCE for having offered the opportunity and the spot 
its journal for the publication his paper. 

“February, 1980, Sam Chan (Proc. Paper 15211). 


Coordinator Engrg., Div. Sci., Louisiana State Univ., P.O. Box 1129, Eunice, 
La. 70535. 
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The writers thank Gupta for his contributions their study. 

The discusser’s comments the maximum frequency content wind records 
essentially agree qualitatively with the writers’ study although somewhat 
severe, the writers’ opinion. 

The use longer rise times suggested the discussion would ideal 
but would add inordinately the computational expense. The discusser’s dynamic 
load factor computations may not recognize the initiation the time history 
analysis from the mean pressure which partially compensates for 
the short rise times. 


Finally, the use initial damping seems sufficiently conservative. 


b 
Errects CONCRETE SHELLS 


Errata 


The following correction should made the original paper: 


“March, 1980, Prodyot Basu and Phillip Gould (Proc. Paper 15259). 

Prof., Dept. Civ. Engrg., Washington Univ., St. Louis, Mo. 63130. 

and Chmn. Civ. Engrg., Washington Univ., Campus Box 1130, St. Louis, 
Mo. 63130. 


April, 1980, Wai-Fah Chen, Hiroyuki Suzuki, and Tse-Yung Chang (Proc. Paper 
15316). 
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River 


Hans and George 


The writers have reviewed the comments provided Elkow and 
have the following observations. 

are agreement with the discusser that the high tensile residual stress 
from welding and the joint geometry are major contributors crack extension, 
particularly when large lack fusion region from partially made weld exists. 
However, are not able ascertain the reasons for his concluding that 
made groove weld would experience cracking, had been erected without 
flaw and then been subjected service temperature drop 

The intersection sound groove welds continuous fillet welds were 
used attach the longitudinal stiffener the girder web common practice 
and also occurred elsewhere the structure question. adverse behavior 
has been experienced such locations when sound groove welds have been 
produced. 

sound groove weld without flaw will not crack result shrinkage 
from longitudinal fillet groove welded connection crossing it, unless poor 
welding practice utilized, e.g., introducing hydrogen the weld. 

The writers are not sure the analytical basis the statements made. 
would appear that the discusser suggesting alternate method analysis. 
However, lack information prevents more detailed discussion and 
evaluation. 

should also noted that the thickness the longitudinal stiffener and 
web plate results residual tensile stresses through the web thickness the 
stiffener-web connection. The longitudinal welds were made the submerged 
arc process and the penetration the weld into the stiffener and web plates 
sufficient cause tensile yield stress through the web thickness the 
longitudinal stiffener connection. 

This condition not severe the web-to-flange connection, because the 
flange plate much thicker than the web plate. This results different 
residual stress distribution was shown Fig. 16. The residual stress distribu- 
tions have impact the fracture behavior the two welded connections 
that were examined the paper. The residual stress gradient the flange 
was factor arresting the running crack that developed the girder web. 

1980, John Fisher, Alan Pense, Hans Hausammann, and George 
Irwin (Proc. Paper 15343). 


Civ. Engrg. and Assoc. Dir., Fritz Engrg. Lab., Building 13, Lehigh Univ., 
Bethlehem, Pa. 18015. 


and Chmn., Metallurgy and Materials Engrg., Lehigh Univ., Bethlehem, Pa. 
18015. 
Asst., Fritz Engrg. Lab., Lehigh Univ., Bethlehem, Pa. 18015. 


Adjunct Prof. Mechanics, Fritz Engrg. Lab., Lehigh University, Bethlehem, Pa. 
18015. 
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Discussion Robert Martin’ and Robert 


The author addresses very practical consideration, namely, the cost 
conducting dynamic seismic analyses which frequently are conducted conjunc- 
tion with other stress calculations. The computer models illustrated, however, 
appear restricted beam elements. Where components cannot modeled 
beams, the author acknowledges that detailed finite element analysis may 
required obtain equivalent beam properties. mention made the 
cost these subsidiary calculations, although the author claims that general 
this economic approach. other instances, the adequacy the models 
illustrated seems rest heavily the fact that the frequency predicted 
anaiysis conservative with respect some known value (perhaps established 
test?) for the physical article. many instances, the analyst does not have 
the advantage such prior knowledge for assessing the adequacy his computer 
model and, therefore, must resort more detailed model, while the same 
time living within the constraints the inevitable budgetary limitations. 

The writers have considerable experience modeling equipment varying 
complexity (some for nuclear applications) requiring dynamic seismic/stress 
analysis. general purpose structural analysis program, e.g., NASTRAN, 
ANSYS, having fairly extensive library finite elements (beam, quad plate, 
spring, concentrated mass, etc.) has not only proven useful, but with 
experience, rather economical. consequence, has been possible 
synthesize, much the same way the author’s paper suggests, basic library 
computer models for various types fans, cabinets, etc. which can 
used repeatedly, thus resulting cost effectiveness. Moreover, single model 
adequate for both dynamic and static analysis since condensation technique, 
Guyan reduction, available facilitating the resonant frequency computations, 
while re-expansion the modes available for the dynamic seismic stress 
portion the analysis. Thus, one model suffices for the entire problem, 
distinction two models the author’s approach. Certain post-processing 
features make possible and economical combine dynamic stress results 
with static stresses arising from dead weight and operating loads. 

1980, McDonald (Proc. Paper 15571). 


Computer Operations, McMahon Engrg. Co., 2322 Bowen Rd., P.O. Box 
Elma, N.Y. 14059. 


Analyst, McMahon Engrg. Co., 2322 Bowen Rd., P.O. Box Elma, N.Y. 
14059. 
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